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Preface 


RANDOMIZED BY FOOLISHNESS 


This is a book that has been written very rapidly in order to respond 
to a perceived need for clarity in the quantitative world. The second 
section of this preface, “How I learned to stop worrying and love 
the CDOs”, is obviously ironic in referring to the hysteria that has 
often characterized accounts of modelling and mathematical finance 
in part of the press and the media, and the demonization of part of the 
market products related to the crisis, such as CDOs and derivatives 
more generally (Brigo et al., 2009b). Let us be clear and avoid any 
misunderstanding: the crisis is very real; it has caused suffering to 
many individuals, families and companies. However, it does not help 
to look for a scapegoat without looking at the whole picture with a 
critical eye. Accounts that try to convince the public that the crisis 
is due mainly to modelling and to the trading of sophisticated and 
obscure products are necessarily partial, and this book is devoted to 
rectifying this perceptive bias. 

Indeed, public opinion has been bombarded with so many clichés 
on derivatives, modelling and quantitative analysis that we feel that 
a book offering a little clarity is needed. And while we are aware 
that this sounds a little Don Quixotesque, we hope the book will 
help to change the situation. In trying to do so, we need to balance 
carefully the perspectives of different readerships. We would like our 
book to be attractive to the relatively general industry and academic 
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public without disappointing the scientific and technically minded 
specialists, and at the same time we do not want our book to be a best- 
seller type of publication full of bashing and negatively provocative 
ideas with very little actual technical content. All of this, while 
keeping windmills at large.! Hence we will be walking the razor’s 
hedge in trying to maintain a balance between the popular account 
and the scientific discourse. 

We are not alone in our attempt to bring clarity.” This book, how- 
ever, does so in an extensive and technical way, showing past and 
present research that is quite relevant in disproving a number of mis- 
conceptions on the role of mathematics and quantitative analysis in 
relation to the crisis. This book takes an extensive technical path, 
starting with static copulas and ending with dynamic loss models. 
Even though our book is short, we follow a long path for credit 
derivatives and multi-name credit derivatives in particular, focusing 
on Collateralized Debt Obligations (CDOs). What are CDOs? To 
describe the simplest possible CDO, say a synthetic CDO on the 
corporate market, we can proceed as follows. 

Weare given a portfolio of, say, 125 names. The names may default 
thus generating losses to investors exposed to those names. In a CDO 
tranche there are two parties: a protection buyer and a protection 
seller. A tranche is a portion of the loss of the portfolio between two 
percentages. For example, the 3—6% tranche focuses on the losses 
between 3% (attachment point) and 6% (detachment point). Roughly 
speaking, the protection seller agrees to pay to the buyer all notional 
default losses (minus the recoveries) in the portfolio whenever they 
occur due to one or more defaults of the entities, within 3% and 6% 
of the total pool loss. In exchange for this, the buyer pays the seller a 
periodic fee on the notional given by the portion of the tranche that 
is still “alive” in each relevant period. 


! The Ingenious Hidalgo Don Quixote of La Mancha, Miguel de Cervantes, 1605 and 1615. 

2 See primarily Shreve (2008) and also, for example, Donnelly and Embrechts (2009). Szego 
(2009) gives a much broader overview of the crisis, with some critical insight that is helpful in 
clarifying the more common misconceptions. 
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In a sense, CDOs look like contracts selling (or buying) insurance 
on portions of the loss of a portfolio. The valuation problem is trying 
to determine the fair price of this insurance. 

The crucial observation here is that “tranching” is a non-linear 
operation. When computing the price (mark to market) of a tranche 
at a point in time, one has to take the expectation of the future tranche 
losses under the pricing measure. Since the tranche is a non-linear 
function of the loss, the expectation will depend on all moments of the 
loss and not just on the expected loss. If we look at the single names 
in the portfolio, the loss distribution of the portfolio is characterized 
by the marginal distributions of the single-name defaults and by the 
dependency among the defaults of different names. Dependence is 
commonly called, with an abuse of language, “correlation”. This is 
an abuse of language because correlation is a complete description 
of dependence for jointly Gaussian random variables, although more 
generally it is not. In general, the term “dependence” should be used 
in lieu of “correlation”. The complete description is either the whole 
multivariate distribution or the so-called “copula function” — that 
is, the multivariate distribution once the marginal distributions have 
been standardized to uniform distributions. 

The dependence of the tranche on “correlation” is crucial. The mar- 
ket assumes that there is a Gaussian Copula connecting the defaults 
of the 125 names. This copula is parameterized by a matrix with 7750 
entries of pairwise correlation parameters. However, when looking 
at a tranche, these 7750 parameters are all assumed to be equal to 
each other. So one has a unique parameter. This is such a drastic 
simplification that we need to make sure it is noticed: 


7750 parameters —> 1 parameter. 


One then chooses the tranches that are liquid on the market for 
standardized portfolios, for which the market price is known as these 
tranches are quoted. The unique correlation parameter is then reverse- 
engineered to reproduce the price of the liquid tranche under exami- 
nation. This is called implied correlation, and once obtained it is used 
to value related products. The problem is that whenever the tranche 
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is changed, this implied correlation also changes. Therefore, if at a 
given time the 3—6% tranche for a 5-year maturity has a given im- 
plied correlation, the 6-9% tranche for the same maturity will have 
a different one. It follows that the two tranches on the same pool 
are priced with two models having different and inconsistent loss 
distributions, corresponding to the two different correlation values 
that have been implied. 

This may sound negative, but as a matter of fact the situation is even 
worse. We will explain in detail that there are two possible implied 
correlation paradigms: compound correlation and base correlation. 
The second correlation is the one that is prevailing in the market. 
However, the base correlation is inconsistent even at a single tranche 
level, in that it prices the 3—6% tranche by decomposing it into the 
0-3% tranche and 0—6% tranche and, for that tranche level, using 
two different correlations (and hence distributions). Therefore, the 
base correlation is already inconsistent at the single tranche level, 
and this inconsistency shows up occasionally in negative losses (i.e. 
in defaulted names resurrecting). 

This is admittedly enough to spark a debate. Even before modelling 
enters the picture, some famous market protagonists have labelled the 
objects of modelling (i.e. derivatives) as being responsible for a lot of 
problems. Warren Buffett (2003), in a very interesting report, wrote: 


[... ] Charlie and I are of one mind in how we feel about derivatives 
and the trading activities that go with them: We view them as time 
bombs, both for the parties that deal in them and the economic system. 
[...] In our view [...] derivatives are financial weapons of mass 
destruction, carrying dangers that, while now latent, are potentially 
lethal. [... ] The range of derivatives contracts is limited only by the 
imagination of man (or sometimes, so it seems, madmen). 


While, when hearing about products such as Constant Proportion 
Debt Obligations (CPDOs) or CDO squared, one may sympathize 
with Mr Buffett, this overgeneralization might be a little excessive. 
Derivatives, when used properly, can be quite useful. For exam- 
ple, swap contracts on several asset classes (interest rates, foreign 
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exchange, oil and other commodities) and related options allow enti- 
ties to trade risks and buy protection against adverse market moves. 
Without derivatives, companies could not protect themselves against 
adverse future movements in the prices of oil, exchange rates, interest 
rates, etc. This is not to say that derivatives cannot be abused. They 
certainly can, and we invite the interested-readers to reason on the 
case of CPDOs? as an example, and to read the whole report by Mr 
Buffett. 

When moving beyond the products and entering the modelling 
issues, one may still find popular accounts resorting to quite colourful 
expressions such as “the formula that killed Wall Street’. Indeed, if 
one looks at popular accounts such as Salmon (2009) or Jones (2009), 
just to take two examples, one may end up with the impression that the 
quantitative finance (“quant”) community has been incredibly naive 
in accepting the Gaussian Copula and implied correlation without 
questioning it, possibly leading to what Warren Buffett calls “mark 
to myth” in his above-mentioned report, especially when applying 
the calibrated correlation to other non-quoted “bespoke tranches”. In 
fact both articles have been written on the Gaussian Copula —a static 
model that is little more than a static multivariate distribution which is 
used in credit derivatives (and in particular CDOs) valuation and risk 
management. Can this simple static model have fooled everyone into 
believing it was an accurate representation of a quite dynamic reality, 
and further cause the downfall of Wall Street banks? While Salmon 
(2009) correctly reports that some of the deficiencies of the model 
have been known, Jones (2009) in the Financial Times wonders why 
no one seemed to have noticed the model’s weaknesses. The crisis 
is considered to have been heavily affected by mathematical models, 
with the accent on “mathematical”. 

This is in line with more general criticism of anything quantitative 
appearing in the news. As an example, the news article “McCormick 


3 See, for example, Torresetti and Pallavicini (2007) and the Fitch Ratings report “First Gener- 
ation CPDO: Case Study on Performance and Ratings’, published before the crisis in April 2007, 
stating “[...] Fitch is of the opinion that the past 10 years by no means marked a high investment 
grade stress in the range of AAA or AA.” 
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Bad Dollars Derive From Deficits Model Beating Quants”, which is 
more focused on the currency markets, informs us that “[...] focus on 
the economic reasons for currency moves is gaining more traction af- 
ter years when traders and investors relied on mathematical models of 
quantitative analysis.” Then it continues with “These tools worked 
during times of global growth and declining volatility earlier this 
decade, yet failed to signal danger before the financial crisis sparked 
the biggest currency swings in more than 15 years. McCormick, using 
macroeconomic and quantitative analyses, detected growing stresses 
in the global economy before the meltdown.” The reader, by looking 
at this, may understand that on one side “mathematical models of 
quantitative analysis” (a sentence that sounds quite redundant) fail in 
times of crises, whereas “macroeconomic and quantitative analyses” 
helped to predict some aspects of the crisis. One is left to ponder 
the different uses of “quantitative” between “mathematical models 
of quantitative analysis” and “macroeconomic and quantitative anal- 
yses”. It is as if mathematics had all of a sudden become a bad word. 
Of course the article aims at saying that macroeconomic analysis 
and fundamentals are of increasing importance and should be taken 
more into account, although in our opinion it does not clearly distin- 
guish valuation from prediction, but some of the sentences used to 
highlight this idea are quite symptomatic of the attitude we described 
above towards modelling and mathematics. 

Another article that brings mathematics and mathematicians 
(provided that is what one means by “math wizards”) into the “blam- 
ing” picture is Lohr (2009) — “Wall Street’s Math Wizards Forgot 
a Few Variables” — which appeared in the New York Times of 12 
September. Also, Turner (2009) has a section entitled “Misplaced 
reliance on sophisticated maths”. 

This overall hostility and blaming attitude towards mathematics 
and mathematicians, whether in the industry or in academia, is the 
reason why we feel it is important to point out the following: the 
notion that even more mathematically oriented quants have not been 


4 By Oliver Biggadike, 24 November 2009, Bloomberg. 
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aware of the Gaussian Copula model’s limitations is simply false, 
as we are going to show, and you may quote us on this. The quant 
and academic communities have produced and witnessed a large 
body of research questioning the copula assumption. This is well 
documented: there is even a book based on a one-day conference 
hosted by Merrill Lynch in London in 2006, well before the crisis, and 
called Credit Correlation: Life After Copulas (Lipton and Rennie, 
2007). This conference had been organized by practitioners. The 
Life After Copulas book contains several attempts to go beyond the 
Gaussian Copula and implied correlation, most of which come from 
practitioners (and a few from academics). But that book is only the 
tip of the iceberg. There are several publications that appeared pre- 
crisis that questioned the Gaussian Copula and implied correlation. 
For example, we warned against the dangers implicit in the use of 
implied correlation in the report “Implied Correlation: A paradigm 
to be handled with care”, that we posted in SSRN in 2006, again well 
before the crisis (Torresetti et al., 2006b). 

Still, it seems that this is little appreciated by some market partic- 
ipants, commentators, journalists, critics, politicians and academics. 
There are still a number of people who think that a formula killed 
Wall Street (see also the discussion in Embrechts, 2009). 

This book brings a little clarity by telling a true story of pre-crisis 
warnings and also of pre-crisis attempts to remedy the drawbacks of 
implied correlation. We do not document the whole body of research 
that has addressed the limits of base correlation and of the Gaussian 
Copula, but rather take a particular path inside this body, based on 
our past research, which we also update to see what our models tell 
us in-crisis. 

To explain our book in a nutshell, we can say that it starts from 
the payoffs of CDOs, explaining how to write them and how they 
work. We then move to the introduction of the inconsistent Gaussian 
Copula model and the related implied correlations, both compound 
and base, moving then to the GPL model: an arbitrage-free dynamic 
loss model capable of consistently calibrating all the tranches across 
attachments and detachments for all the maturities at the same time. 
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En passant, we also illustrate the Implied Copula, a method that 
can consistently account for CDOs with different attachment and 
detachment points but not for different maturities, and the Expected 
Tranche Loss (ETL) surface, a model-independent approach to CDO 
price interpolation. 

We will see that, already pre-crisis, both the Implied Copula and 
the dynamic loss model imply modes far down the right tail of the 
loss distribution. This means that there are default probability clusters 
corresponding to joint default of a large number of entities (sectors) 
of the economy. 

The discussion is abundantly supported by market examples 
throughout history. We cannot stress enough that the dangers and 
critics we present to the use of the Gaussian Copula and of implied 
correlation, and the modes in the tail of the loss distribution obtained 
with consistent models, had all been published by us, among others, 
in 2006, well before the crisis. 

Despite these warnings, the Gaussian Copula model is still used in 
its base correlation formulation, although under some possible ex- 
tensions such as random recovery. The reasons for this are complex. 
First, the difficulty of all the loss models, improving the consistency 
issues, in accounting for single-name data and to allow for single- 
name sensitivities. This is due to the fact that if we model the loss of 
the pool directly as an aggregate object, without taking into account 
single defaults, then the model sensitivities to single-name credit in- 
formation are not in the picture. In other terms, while the aggregate 
loss is modelled so as to calibrate satisfactorily indices and tranches, 
the model does not see the single-name defaults but just the loss dy- 
namics as an aggregate object. Therefore partial hedges with respect 
to single names are not possible. As these issues are crucial in many 
situations, the market practice remains with base correlation. Further- 
more, even the few models achieving single-name consistency have 
not been developed and tested enough to become operational on a 
trading floor or on a large risk management platform. Indeed, a fully 
operational model with realistic run times and numerical stability 
across a large range of possible market inputs would be more than 
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a prototype with some satisfactory properties that has been run in 
some “‘off-line” studies. Also, when one model has been coded in the 
libraries of a financial institution, changing the model implies a long 
path involving a number of issues that have little to do with modelling 
and more to do with IT problems, integration with other systems, etc. 
Therefore, unless a new model appears to be really promising and 
extremely convincing in all its aspects, there is reluctance to adopt it 
on the trading floor or for risk management systems. 

Overall we conclude that the modelling effort in this area of the 
derivatives market is unfinished, partly for the lack of an operationally 
attractive single-name consistent dynamic loss model, and partly 
because of the diminished investment in this research area, but the 
fact that the modelling effort is unfinished does not mean that the 
quant community has been unaware of model limitations, as we 
abundantly document, and, although our narrative ends with an open 
finale, we still think it is an entertaining true story. 


HOW I LEARNED TO STOP WORRYING 
AND LOVE THE CDOs 


We cannot close this preface without going back to the large picture, 
and ask the more general question: Is the crisis due to poor modelling? 

As we have seen, the market has been using simplistic approaches 
for credit derivatives, but it has also been trying to move beyond those 
approaches. However, we should also mention that CDOs are divided 
into two categories: cash and synthetics. Cash CDOs involve hun- 
dreds or even thousands of names and have complex path-dependent 
payouts (“waterfalls”). Even so, cash CDOs are typically valued by 
resorting to single homogeneous default-rate scenarios or very prim- 
itive assumptions, and very little research and literature is available 
on them. Hence, these are complex products with sophisticated and 
path-dependent payouts that are often valued with extremely sim- 
plistic models. Synthetic CDOs are the ones that we described in this 
preface and will be addressed in this book. They have more simple 
and standardized payouts than the cash CDOs but are typically valued 
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with more sophisticated models, given the larger standardization and 
the ease in finding market quotes for their prices. Synthetic CDOs 
on corporates are epitomized by the quoted tranches of the stan- 
dard pools DJ-iTraxx (Europe) and CDX (USA). However, CDOs, 
especially cash CDOs, are available on other asset classes, such as 
loans (CLO), Mortgage Portfolios, leading to Residential Mortgage 
Backed Portfolios (RMBS) and Commercial Mortgage Backed Se- 
curities (CMBS), and on and on. For many of these CDOs, and 
especially RMBSs which are quite related to the asset class that trig- 
gered the crisis, the problem is in the data rather than in the models. 
Bespoke corporate pools have no data from which to infer default 
“correlation” and dubious mapping methods are used. At times, data 
for valuation in mortgages CDOs (RMBSs and CDOs of RMBSs) 
are dubious and can be distorted by fraud.’ 

At times it is not even clear what is in the portfolio: the authors have 
seen offering circulars of a RMBS on a huge portfolio of residential 
mortgages where more than the 2% of properties in the portfolio 
were declared to be of unknown type. What inputs can we give to 
the models if we do not even know the kind of residential property 
underlying the derivative product? 

All this is before modelling. Models obey a simple rule that is pop- 
ularly summarized by the acronym GIGO (Garbage In — Garbage 
Out). As Charles Babbage (1791-1871) famously put it: 


On two occasions I have been asked [by members of Parliament], 
“Pray, Mr. Babbage, if you put into the machine wrong figures, will 
the right answers come out?” I am not able rightly to apprehend the 
kind of confusion of ideas that could provoke such a question. 


So, in the end, is the crisis due to the inadequacy of the models? 
Is the crisis due to the pride of quantitative analysts and academics 
and an unawareness of the limitations of the models? 

We show in this book that quants have been aware of the limitations 
and of extreme risks before the crisis. Lack of data or fraud-corrupted 


5 See, for example, the FBI Mortgage Fraud Report, 2007, www.fbi.gov/publications/fraud/ 
mortgage_fraud07.htm. 
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data, the fragility in the “originate to distribute” system, liquidity 
and reserves policies, regulators’ lack of uniformity, excessive lever- 
age and concentration in real estate investment, poor liquidity risk 
management techniques, accounting rules and excessive reliance on 
credit rating agencies are often factors not to be underestimated. This 
crisis is a quite complex event that defies witch-hunts, folklore and 
superstition. Methodology certainly needs to be improved but blam- 
ing just the models for the crisis appears, in our opinion, to be the 
result of a very limited point of view. 


London, Milan, Madrid, Pavia and Venice, 1 December 2009. 


Damiano Brigo, Andrea Pallavicini and Roberto Torresetti. 
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| Notation and List of Symbols 


Probabilities and expectations 


e Q: Risk-neutral probability. 

e E: Risk-neutral expectation. 

e px: Risk-neutral probability density function of the random vari- 
able X: 


d 
px(u) := mL <u} 
u 


e ox: Risk-neutral characteristic function of the random variable X, 
namely the Fourier transform of the probability density function: 


gx (u) := Eo[ e"* ] = F[ px 1 (u) 


e hy: Risk-neutral intensity of the process X, namely the density of 
the absolutely continuous compensator of X: 
El Xusta] — Xu 


t 
; d 
Ax(t):= lim f De du, hyx(t):= T Ax(t) 


e y: Probability density function of the standard n-variate Gaussian 
distribution with correlation matrix C: 


1 LTC- 
a(x, C) = —3x' Cx 
n(x, C) Onn’ 


In particular for the case n = 1, 2 we define: 


1 
P(x) = p, 1), pax, X2, €) = Q2 (on, x2), e J) 


xxvi Notation and List of Symbols 


e : Cumulative distribution function of the standard n-variate Gaus- 
sian distribution with correlation matrix C: 


P, (Xx, o= f ai Only, C) dy, ...dyn 


In particular, for the case n = 1, 2 we define: 


1 
PO) Oily Di -Po E (1.20) (è D 


e f: Probability density function of the Student t-distribution with 
d degrees of freedom: 


j 
f. ( ) r (5(d + 1)) (1 ae ae 
x) i= Se 
i Vax I (4d) 
e F: Cumulative distribution function of the Student t-distribution 
with d degrees of freedom: 


d 


F(x) := J Jfa) dy 


Default events and Poisson processes 


e t: Default time of one name within the pool. 
e p: Default probability of one name within the pool: 


pit) = Qiu <t} 


e A: Possibly stochastic default intensity (and then hazard rate). 
e A: Possibly stochastic cumulated default intensity (and then hazard 
function): 


A(t) =) Au) du 
0 


e N: Poisson process with a possibly stochastic intensity À. 
e 1: Indicator function: 


; 1l,t>T 
Mon = Jo per 
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Counting and loss processes 


e M: Number of elements of a pool of names. For the main DJ-iTraxx 
and CDX indices, M is equal to 125. 

e A: Tranche attachment point. 

B: Tranche detachment point. 

C: Pool’s counting process: 


M 
C, := 2 lirst) 
i=l 


e Č: Pool’s normalized counting process: 


= 1 
C; = ie 


e R: Recovery rate of one name within the pool. 
e Lo: Loss given default one name within the pool: 


Leo; = 1- R; 


e L: Pool’s loss process: 


M 


L; = 5 lir <1 Loo; 


i=l 
e L: Pool’s normalized loss process: 


£ 1 
L, = aa 
e L4-8: Pool’s normalized loss process of a CDO tranche with at- 


tachment point A and detachment point B: 


1 
B-A 


Lee [Li - A) acta) + (B - Aliza] 


e £4:8-: Pool’s normalized loss of a CDO tranche with attachment 
point A and detachment point B under the large pool assumptions. 
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e f(t, A, B): Pool’s normalized expected loss of a CDO tranche with 
attachment point A and detachment point B: 


f(t, A, B) := Eo| £4? | 


e g(t, B): Pool’s normalized expected loss of a CDO equity tranche 
with detachment point B. 


g(t, B) := B f (t,0, B) 


Copula models 


e p”: Base correlation of a CDO equity tranche with detachment 
point B. 

e a g: Compound correlation of a CDO tranche with attachment 
point A and detachment point B. 

e S: Systemic factor affecting default times of all names in one- 

factor copula models. In the Gaussian Copula model it is a Gaussian 

random variable, while in the Double-t Copula model it is a Student 

t-variable. 

Y: Idiosyncratic factor affecting just one name in one-factor copula 


models. In the Gaussian Copula model it is Gaussian random vari- 

able, while in the Double-t Copula model it is a Student t-variable. 
e X: Random variable entering the definition of copula models. For 

example, in the Gaussian copula, for each name i we get: 


X; := 07! (1 = e^) 


If we consider a one-factor version of the model, then the following 
equation relates the systemic and idiosyncratic factors to the copula 
with p; as correlation parameter for the ith name of the pool: 


Xi =4pPi S +vy1-— pi Y; 
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Dynamic loss models 


e M’: Inverse of the minimum jump for the pool’s loss process. The 
following constraints hold: 


M>M>0 
e R: Average recovery rate: 
R:=1- Fol Li | 
Eo[ C, | 


e a: Number of defaults triggered by the jump of the jth Poisson 
processes contributing to the GPL model. 

e a: Fraction of defaults triggered by the jump of the jth Poisson 
processes contributing to the GPL model: 


sah 
aj := We 


e Z: GPL driving process: 


Z= X ajNj® 


j=l 
where N; is a Poisson process with a possibly stochastic intensity 
Àj. 
e oz: Characteristic function of the GPL driving process: 


gz (u) = Eo| exp | > Aj(t)(e""” — 1) 
j=l 


e ĉ: The first time when the GPL driving process reaches or exceeds 
the number of names of the pool: 


îĉ = inf {t : Z, > M} 


e Y; A deterministic non-decreasing function defining in the GPL 
model the pool’s loss process: 


Y : NU {0} —> [0,1], L; := Vz(Z,) 
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Product’s payoffs and prices 


e D(t, T): Discount factor at time t for maturity T. 
e P(t, T): Default-free zero-coupon bond price at time ¢ for maturity 
T: 


P(t, T):=E,[ Dd, T)] 


If interest rates are taken as deterministic, then P(t, T) = D(t, T). 
e {7|,..., Tp = T}: Payment dates of the periodic fee (spread) S 
payed by the premium leg of a CDO. 
e 5: Year fraction appearing in the CDO premium leg. 
e S: Running spread payed by the premium leg of a CDO index: 


ARA Dit, dL 
E| Dh. 8DE, T — Čr) | 


Si = 


e DixLec: Default leg payoff of a CDO index starting at time t with 
maturity 7: 


T 
DerLeo(t) := J D(t,u)dL, 
t 
e PremumLec: Premium leg payoff of a CDO index starting at time t 


with maturity T: 


b 
PrewusLes(t) = S, XC 8;D(t, TOA — Čr) 
i=l 


e Nov: Net present value of a CDO index starting at time ¢ with 
maturity 7: 


Nov(t) = EL PremumLec(t) = DerLec(t) ] 


e U4-8: Upfront premium payed to enter a CDO tranche with at- 
tachment point A and detachment point B. 
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e S^?: Running spread payed by the premium leg of a CDO tranche: 


B| f De, wale? | — uh? 
E,[ Era DG, TU — ZF”) | 


e DirLec, g: Default leg payoff of a CDO tranche with attachment 
point A and detachment point B starting at time t with maturity T: 


A,B ._ 
SAB = 


T 
DerLec 4. p(t) =| D(t, u)dLA? 
t 


e Dvı4 g: Annuity of a CDO tranche with attachment point A and 
detachment point B starting at time t with maturity T: 


b 
Duwisa(t) := X 5: DE, TU — Ly”) 
i=1 
o PremumLec4,g: Premium leg payoff of a CDO tranche with attach- 
ment point A and detachment point B starting at time t with ma- 
turity T: 
PremumLec 4 p(t) = pee + S&P Dvo p(t) 


e N»v4,g: Net present value of a CDO tranche with attachment point 
A and detachment point B starting at time t with maturity T: 


Nova, g(t) = E; [ PrewumLec4, g(t) = DeLes4, g(t) | 
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Introduction: 
Credit Modelling Pre- and 
In-Crisis 


This book aims to show the limits of popular models or pseudo- 
models (mostly quoting mechanisms with modelling semblance) 
that, in recent years, have been extensively used to mark to mar- 
ket and risk manage multi-name credit derivatives. We present a 
compendium of results that we first published in 2006 before the cri- 
sis pointing out the dangers in the modelling paradigms used at the 
time in the market, and showing how the situation has even worsened 
subsequently by analysing more recent data. We also point out that 
the current paradigm had been heavily criticized before the crisis, 
referring to our work and the works of other authors addressing the 
main limitations of the current market paradigm well before popu- 
lar accounts such as Salmon (2009) appeared. The problems of the 
current paradigm include: 


e An unrealistic Gaussian Copula assumption and the flattening of 
7750 pairwise dependence parameters into one. 

e Lack of consistency of the implied correlation market models with 
more than one tranche quote at the time. 

e Occasional impossibility of calibration even of single tranches, or 
possibility to obtain negative expected tranched losses violating 
the arbitrage-free constraints. 

e Lack of an implied loss distribution consistent with market CDO 
tranche quotes for a single maturity. 

e Lack of a loss distribution dynamics consistent with CDO tranche 
quotes on several maturities. 


2 Credit Models and the Crisis 


e Lack of credit spread volatility, resulting for example in heavy 
consequences on the valuation of counterparty risk on CDS under 
wrong way risk. 


In this respect we will introduce examples of models published before 
the crisis that partly remedy the above deficiencies. All the discussion 
is supported by examples based on market data, pre- and in-crisis. 
In addressing these issues we adopt the following path through the 
different methodologies. 


1.1 BOTTOM-UP MODELS 


A common way to introduce dependence in credit derivatives mod- 
elling is by means of copula functions. A typically Gaussian Copula 
is postulated on the exponential random variables triggering defaults 
of the pool names according to the first jumps of Poisson processes. In 
general, if one tries to model dependence by specifying dependence 
across single default times, one is in a so-called “bottom-up” frame- 
work, and the copula approach is typically within this framework. 
Such a procedure cannot easily be extended to a fully dynamical 
model in general. We cannot do justice to the huge copula literature 
in credit derivatives here; we only mention that there have been at- 
tempts to go beyond the Gaussian Copula introduced in the CDO 
world by Li (2000) and leading to the implied (base and compound) 
correlation framework, some important limits of which have been 
pointed out in Torresetti et al. (2006b). Li and Hong Liang (2005) 
also proposed a mixture approach in connection with CDO squared. 
For the results on sensitivities computed with the Gaussian Copula 
models, see for example Meng and Sengupta (2008). 

An alternative to copulas in the bottom-up context is to insert de- 
pendence among the default intensities of single names — see, for ex- 
ample, the paper by Chapovsky et al. (2007). Joshi and Stacey (2006) 
resort to modelling business time to create default correlation in oth- 
erwise independent single-name defaults, resorting to an “intensity 
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gamma” framework. Similarly, but in a firm value-inspired context, 
Baxter (2007) introduces Levy firm value processes in a bottom- 
up framework for CDO calibration. Lopatin (2008) introduces a 
bottom-up framework that is also effective in the CDO context, hav- 
ing single-name default intensities that are deterministic functions of 
time and of the pool default counting process, then focusing on hedge 
ratios and analysing the framework from a numerical performances 
point of view, showing this model to be interesting even if lacking 
explicit modelling of single-name credit spread volatilities. 

It is worth noticing that copula models are usually implemented 
with deterministic credit spreads. Credit spread volatility is assumed 
to be zero even if both historical (Hull and White, 2003) and implied 
(Brigo, 2005, 2006) CDS volatilities attain values above 50%. In 
first-to-default baskets, which are similar in a way to equity-tranches 
of CDOs (see, for example, Brigo and Mercurio, 2006), spread dis- 
persion combined with lack of credit spread volatility may result in 
unintuitive features of the model prices with respect to the copula cor- 
relation parameter. This has been highlighted in the different context 
of counterparty risk for credit default swaps by Brigo and Chourdakis 
(2009), who show that credit spread volatility can be quite relevant, 
to the point that neglecting it changes the wrong way risk profile 
of CDS pricing under counterparty risk. This is further highlighted 
by Brigo and Capponi (2008), while Brigo and Pallavicini (2007, 
2008), Brigo and Bakkar (2009) and Brigo et al. (2009) model credit 
spread volatility explicitly, showing that credit spread volatility has 
a relevant impact also on counterparty risk for other asset classes, in- 
cluding commodities and interest rate products. The stochastic credit 
spread models used in these works build on the work of Brigo and 
Alfonsi (2005), Brigo and Cousot (2006) and Brigo and El-Bachir 
(2008) on CDS options. 

Going back to bottom-up models in the context of CDOs, Albanese 
et al. (2006) introduce a bottom-up approach based on structural 
model ideas that can not only be made consistent with several inputs 
under historical and pricing measures, but also manages to calibrate 
CDO tranches. 


4 Credit Models and the Crisis 


1.22 COMPOUND CORRELATION 


Building on Torresetti et al. (2006b), in the context of bottom-up 
models, we start with the net present value (NPV) of synthetic Col- 
lateralized Debt Obligations (CDOs) tranches on pools of corporate 
credit references in its original layout: the compound correlation 
framework. 

We highlight two of the major weaknesses of the compound cor- 
relation: 


e Lack of robustness of the compound correlation framework in 
view of the non-invertibility of mainly the 10-year maturity 
DJ-iTraxx 6-9% and CDX 7-10% tranches and, more recently, 
the non-invertibility of mainly the 10-year maturity DJ-iTraxx 12— 
22% and CDX 10-15% tranches. 

e More importantly from a practical standpoint, we highlight the 
typical non-smooth behaviour of the compound correlation and 
the resulting difficulties in pricing bespoke CDO tranches. 


1.3 BASE CORRELATION 


We then introduce the next step the industry took (see, for example, 
McGinty and Ahluwalia, 2004), namely the introduction of base 
correlation, as a solution to both problems, given the fact that: 


e the resulting map is much smoother, thus facilitating the pricing 
of bespoke tranche spreads from liquid index tranches; 

e until early 2008 the heterogeneous pool one-factor Gaussian 
Copula base correlation had been consistently invertible from 
index market tranche spreads. 


Nevertheless we will expose some of the known remaining weak- 
nesses of the base correlation framework: 


(1) Depending on the interpolation technique being used, tranche 
spreads could not be arbitrage free. In fact, for senior tranches it 
may well be that the expected tranche loss plotted versus time is 
initially decreasing. 
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(2) The impossibility of inverting correlation for senior AAA and 
super-senior tranches. This problem arose only recently follow- 
ing the eagerness of market participants to buy protection from 
systemic risk exposure almost at any cost. 

(3) Inconsistency at single tranche valuation level, as two compo- 
nents of the same trade are valued with models having two dif- 
ferent parameter values. 

(4) Last, but not least, flattening information on 7750 pairwise 
correlation parameters into a single one for each equity tranche 
trade. 


As an explanation for weakness (1) we point to the fact re- 
lated to item (3), namely that this arises because the NPV of each 
tranche is obtained by computing the expected tranche loss and 
outstanding notional under two different distributions (the distri- 
bution corresponding to the attachment base correlation and that 
corresponding to the detachment base correlation) so that base 
correlation is already an inconsistent notion at the single tranche 
level. 

As an explanation for weakness (2), we point to the fact that 
the deterministic recovery assumption, while being computationally 
very convenient, does not allow us to capture the more recent mar- 
ket conditions. This has been addressed in the implied correlation 
framework by Amraoui and Hitier (2008) and Krekel (2008). How- 
ever, even with this update, the base correlation remains flawed and 
may still lead to negative loss distributions. 

Base correlation, with updates and variants, remains to this day 
the main pricing method for synthetic corporate CDOs, regardless 
of the body of research criticizing it that we hint at above and 
below. 


1.4 IMPLIED COPULA 


We next summarize the concept of Implied Copula (introduced by 
Hull and White (2006) as the “perfect copula”) as a non-parametric 
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model by which to deduce, from a set of market CDO spreads span- 
ning the entire capital structure, the shape of the risk-neutral pool loss 
distribution. The general use of flexible systemic factors was later 
generalized and vastly improved by Rosen and Saunders (2009), 
who also discuss the dynamic implications of the systemic factor 
framework. Factors and dynamics are also discussed in Inglis et al. 
(2008). 

Our calibration results, based on the Implied Copula — already 
seen in Torresetti et al. (2006c) — point out that a consistent loss 
distribution across tranches for a single maturity features modes in 
the tail of the loss distribution. These probability masses on the 
far right tail imply default possibilities for large clusters (possibly 
sectors) of names of the economy. These results had been published 
originally in 2006 on ssrn.com. We will report such features here and 
will find the same features again by following a completely different 
approach below. 

Here we highlight the persistence of the modes (bumps) in the 
right tail of the implied loss distribution: 


e through time, via historical calibrations; 

e through regions, comparing the results of the historical calibration 
to the DJ-iTraxx and the CDX; and 

e through the term structure comparing the results of the calibration 
to different maturities. 


The Implied Copula can calibrate consistently across the capital 
structure, but not across maturities, as it is a model that is inher- 
ently static. The next step thus consists in introducing a dynamic 
loss model. This moves us into the so-called top-down framework 
(although dynamic approaches are also possible in the bottom-up 
context, as we have seen in some of the above references). But be- 
fore analysing the top-down framework in detail, we will make a 
quick diversion for a model-independent approach to CDO tranches 
pricing and interpolation. 
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15 EXPECTED TRANCHE LOSS SURFACE 


The expected tranche loss (ETL) for different detachment points and 
maturities can be viewed as the basic brick on which the components 
of synthetic CDO formulas are built with linear operations (but under 
some non-linear constraints). We explain in detail how the payoffs of 
credit indices and tranches are valued in terms of expected tranched 
losses (ETLs). This methodology, first illustrated pre-crisis in 
Torresetti et al. (2006a), reminds us of Walker’s earlier work (2006) 
and of the formal analysis of the properties of expected tranche loss 
in connection with no arbitrage in Livesey and Schlégl (2006). 

ETLs are natural quantities that can be implied from market data, 
and no-arbitrage constraints on ETLs as attachment points and ma- 
turities change are briefly introduced. As an alternative to the incon- 
sistent notion of implied correlation illustrated earlier, we consider 
the ETL surface, built directly from market quotes given minimal 
interpolation assumptions. We check that the type of interpolation 
does not interfere excessively with the results. Instruments bid/ask 
spreads enter our analysis, contrary to Walker’s (2006) earlier work 
on the ETL implied surface. By doing so we find less violations of 
the no-arbitrage conditions. 

We also mention some further references that appeared later and 
dealt with evolutions of this technique: Parcell and Wood (2007), 
again pre-crisis, consider carefully the impact of different kinds of 
interpolation, whereas Garcia and Goossens (2007) compare the ETL 
between the Gaussian Copula and Lévy models. 

In general the ETL implied surface can be used to value tranches 
with non-standard attachments and maturities as an alternative to 
implied correlation. However, deriving hedge ratios as well as ex- 
trapolation may prove difficult. Also, the ETL is not really a model 
but rather a model-independent stripping algorithm, although the 
particular choice of interpolation may be viewed as a modelling 
choice. Eventually the ETL is not helpful for pricing more advanced 
derivatives such as tranche options or cancellable tranches, because 
the ETL does not specify an explicit dynamics for the loss of the 
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pool. To that we turn now, by looking at the top-down dynamic loss 
models. 


1.6 TOP (DOWN) FRAMEWORK 


One could completely give up single-name default modelling and 
focus on the pool loss and default-counting processes, thus consid- 
ering a dynamical model at the aggregate loss level, associated with 
the loss itself or to some suitably defined loss rates. This is the “‘top- 
down” approach — see, for example, Bennani (2005, 2006); Giesecke 
and Goldberg (2005); Sch6nbucher (2005); Di Graziano and Rogers 
(2005); Brigo et al. (2006a, 2006b); Errais et al. (2006); Lopatin and 
Misirpashaev (2007); among others. The first joint calibration re- 
sults of a dynamic loss model across indices, tranches attachments 
and maturities (available in Brigo et al., 2006a), show that even 
a relatively simple loss dynamics, like a capped generalized Pois- 
son process, suffices to account for the loss distribution dynamical 
features embedded in market quotes. This work also confirms the 
Implied Copula findings of Torresetti et al. (2006c), showing that 
the loss distribution tail features a structured multi-modal behaviour, 
implying non-negligible default probabilities for large fractions of 
the pool of credit references, showing the potential for high losses 
implied by CDO quotes before the beginning of the crisis. Cont 
and Minca (2008) use a non-parametric algorithm for the calibration 
of top models, constructing a risk-neutral default intensity process 
for the portfolio underlying the CDO, looking for the risk-neutral 
loss process “closest” to a prior loss process using relative entropy 
techniques. See also Cont and Savescu (2008). 

However, in general to justify the “down” in “top-down” one needs 
to show that from the aggregate loss model one can recover a poste- 
riori consistency with single-name default processes when they are 
not modelled explicitly. Errais et al. (2006) advocate the use of ran- 
dom thinning techniques for their approach; see also Halperin and 
Tomecek (2008), who delve into more practical issues related to ran- 
dom thinning of general loss models, and Bielecki et al. (2008) who 
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also build semi-static hedging examples and consider cases where 
the portfolio loss process may not provide sufficient statistics. 

Still, it is not often clear for specific models whether a fully consis- 
tent single-name default formulation is possible, given an aggregate 
model as the starting point. 

There is a special “bottom-up” approach that can lead to distinct 
and rich loss dynamics. This approach is based on the common 
Poisson shock (CPS) framework, reviewed in Lindskog and McNeil 
(2003). This approach allows for more than one defaulting name in 
small time intervals, contrary to some of the above-mentioned “top- 
down” approaches. In the “bottom-up” language, one sees that this 
approach leads to a Marshall-Olkin Copula linking the first jump 
(default) times of single names. In the “top-down” language, this 
model looks very similar to the GPL model in Brigo et al. (2006a) 
when one does not cap the number of defaults. 

The problem of the CPS framework is that it allows for repeated 
defaults, which is clearly wrong as any name could default more than 
once. 

In the credit derivatives literature, the CPS framework has been 
used, for example, in Elouerkhaoui (2006) — see also references 
therein. Balakrishna (2006) introduces a semi-analytical approach 
allowing again for more than one default in small time intervals and 
hints at its relationship with the CPS framework, and also shows 
some interesting calibration results. Balakrishna (2007) then gener- 
alizes this earlier paper to include delayed default dependence and 
contagion. 


1.7 GPL AND GPCL MODELS 


Brigo et al. (2007) address the repeated default issue in CPS by 
controlling the clusters default dynamics to avoid repetition. They 
calibrate the obtained model satisfactorily to CDO quotes across at- 
tachments and maturities, but the combinatorial complexity for a non- 
homogeneous version of the model is forbidding, so that the resulting 
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GPCL approach is hard to use successfully in practice when taking 
single names into account. 

Still, in the context of the present book, the Generalized-Poisson 
Loss (GPL) and the Generalized-Poisson Cluster Loss (GPCL) mod- 
els will be useful to show how a loss distribution dynamics consistent 
with CDO market quotes should evolve. 

In this book we summarize the GPL model, leaving aside the GPCL 
model. As explained above, the GPL is a dynamical model for the 
loss, and is able to reprice all tranches and all maturities at the same 
time. We employ here a variant that models the loss directly rather 
than the default-counting process plus recovery. The loss is modelled 
as the sum of independent Poisson processes, each associated with 
the default of a different number of entities, and capped at the pool 
size to avoid infinite defaults. The intuition of these driving Poisson 
processes is that of defaults of sectors, although the amplitudes of 
the sectors vary in our formulation of the model pre- and in-crisis. 
In the new model implementation in-crisis for this book, we fix the 
amplitude of the loss triggered by each cluster of defaults a priori, 
without calibrating it, as we did in our earlier GPL work. This makes 
the calibration more transparent and the calibrated intensities of the 
default of sectors easier to interpret. We point out, however, that the 
precise default of sectors is made rigorous only in the GPCL. 

We highlight how the GPL model is able to reproduce the tail 
multimodal feature that the Implied Copula proved was indispensable 
for accurately repricing the market spreads of CDO tranches on a 
single maturity. 

We also refer to the later related results of Longstaff and Rajan 
(2008), that point in the same direction but add a principal component 
analysis on a panel of CDS spread changes, with some more com- 
ments on the economic interpretation of the default clusters being 
sectors. 

Incidentally, we draw the reader’s attention to the default history, 
pointing to default clusters being concentrated in a relatively short 
time period (a few months) like the thrifts in the early 1990s at 
the height of the loan and deposit crisis, the airliners after 2001 
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and the autos and financials more recently. In particular, from 
7 September 2008 to 8 October 2008 — a time window of one month — 
we witnessed seven credit events occurring to major financial entities: 
Fannie Mae, Freddie Mac, Lehman Brothers, Washington Mutual, 
Landsbanki, Glitnir and Kaupthing. The Fannie Mae and Freddie 
Mac conservatorships were announced on the same date (7 Septem- 
ber 2008) and the appointment of a “receivership committee” for 
the three Icelandic banks (Landsbanki, Glitnir and Kaupthing) was 
announced between 7 and 8 October. 

Moreover, Standard and Poors issued a request for comments re- 
lated to changes in the rating criteria of corporate CDOs.’ Thus 
far agencies have been adopting a multifactor Gaussian Copula ap- 
proach to simulate the portfolio loss in the objective measure. S&P 
proposed changing the criteria so that tranches rated AAA should be 
able to withstand the default of the largest single industry in the asset 
pool with zero recoveries. We believe that this goes in the direction 
of modelling the loss in the risk-neutral measure via GPL-like pro- 
cesses, given that the proposed changes to S&P’s rating criteria imply 
admitting — as a stressed but plausible scenario — the possibility that a 
cluster of defaults in the objective measure exists. See also Torresetti 
and Pallavicini (2007) for the specific case of Constant Proportion 
Debt Obligations (CPDOs). 

We finally comment more generally on the dynamical aggregate 
models and on their difficulties to lead to single-name hedge ratios 
when trying to avoid complex combinatorics. The framework thus 
remains incomplete to this day, because obtaining jointly tractable 
dynamics and consistent single-name hedges that can be realistically 
applied in a trading floor remains a problem. We provided above some 
references for the latest research in this field. We emphasize, however, 
that even a simple dynamical model like our GPL or the single- 
maturity Implied Copula is enough to appreciate that the market 
quotes were implying the presence of large default clusters with 


' See “Request for Comment: Update to Global Methodologies and Assumptions for Corporate 
Cash Flow CDO and Synthetic CDO Ratings”, 18 March 2009, Standard & Poors. 
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non-negligible probabilities well in advance of the credit crisis, as 
we documented in 2006 and early 2007. 

Finally, it is important to point out that most of the above discus- 
sion and references (with very few exceptions) centre on corporate 
CDOs, and mostly on synthetic CDOs, and little literature is avail- 
able for valuation of CDOs on other asset classes, with possibly 
complex waterfalls and prepayment risk that are cash rather than 
synthetic, including Collateralized Loan Obligations (CLOs), Resi- 
dential Mortgage-Backed Securities (RMBSs) and CDOs of RMBSs, 
which are more related to the asset class that triggered the crisis. 
For many such deals the main problem is often the data. The few 
works on this area include Jaeckel (2008) and Papadopoulos and Tan 
(2007). 


1.8 STRUCTURE OF THE BOOK 


The book is structured as follows: Chapter 2 introduces the market 
quotes we are looking at and provides general discounted payoffs, 
arbitrage-free prices and spread formulas for CDS indices and CDO 
tranches. 

Chapter 3 introduces the Gaussian Copula model, in its different 
formulations concerning homogeneity and finiteness, and then illus- 
trates the notions of implied correlation from CDO tranche quotes. 
The two paradigms of base correlation and compound correlation are 
explained in detail. Existence and uniqueness of implied correlation 
are discussed on a number of market examples, highlighting the pros 
and cons of compound and base correlations, and the limitations in- 
herent in these concepts. In particular, Section 3.9 summarizes issues 
with implied correlations, pointing out the danger for arbitrage when 
negative expected tranche losses appear, and the lack of consistency 
across capital structure and maturity. The first inconsistency is then 
addressed in Chapters 4 and 5, with the Implied Copula, illustrated 
with a number of studies throughout a long period, whereas both 
inconsistencies are addressed in Chapter 6, where our fully-fledged 
GPL dynamic loss model is illustrated pre-crisis. En passant, in 
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Chapter 5 we introduce and analyse the notion of Expected Tranche 
Loss (ETL), a quantity that can be implied in a model-independent 
way from index and CDO tranche quotes and can be used for the 
interpolation of CDO quotes across attachments and maturities. 

All these paradigms are then analysed in-crisis in Chapter 7, while 
the final discussion, including the reasons why implied correlation 
is still used despite all its important shortcomings, are given in 
Chapter 8. In particular, the need for hedge ratios with respect to 
single names, random recovery modelling and speed of calibration 
remain issues that are hard to address jointly outside the base corre- 
lation framework. 


2 
Market Quotes 


For single names our reference products will be credit default swaps 
(CDSs). The most liquid multi-name credit instruments available 
in the market are instead credit indices and CDO tranches (e.g. 
DJ-iTraxx, CDX). We discuss them in the following. 

The procedure for selecting the standardized pool of names is the 
same for the two indices. Every six months a new series is rolled out at 
the end of a polling process managed by Markit, where a selected list 
of dealers contributes the ranking of the most liquid CDSs. All credit 
references that are not investment grade are discarded. Each surviving 
credit reference underlying a CDS is assigned to a sector. Each sector 
is contributing a predetermined number of credit references to the 
final pool of names. The rankings of the various dealers for the 
investment grade names are put together to rank the most liquid 
credit references within each sector. 


2.1 CREDIT INDICES 


The index is given by a pool of names 1, 2,..., M, typically M = 
125, each with notional 1/M so that the total pool notional is 1. 
The index default leg consists of protection payments corresponding 
to the defaulted names of the pool. Each time one or more names 
default, the corresponding loss increment is paid to the protection 
buyer, until final maturity T = T, arrives or until all the names in the 
pool have defaulted. 

In exchange for loss increase payments, a periodic premium with 
rate S is paid from the protection buyer to the protection seller, 
until final maturity Tp. This premium is computed on a notional that 
decreases each time a name in the pool defaults, and decreases of an 
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amount corresponding to the notional of that name (without taking 
out the recovery). 

We denote with L, the portfolio cumulated loss and with C, the 
number of defaulted names up to time ¢ divided by M. Since, at 
each default, part of the defaulted notional is recovered, we have 
0 < L, < Č, < 1. The discounted payoff of the two legs of the index 
is given as follows: 


T 
DerLeo(0) := | D(O, t) dL; 
0 


b 
PrevunLes(0) := So X 6: DO, Ti) — Cn) 
i=1 
where D(s, t) is the discount factor (often assumed to be determin- 
istic) between times s and ¢ and 6; = T; — T;_; is the year fraction. 
In the second equation the actual outstanding notional in each 
period would be an average over [T;—1, T;], but we replaced it with 
1 — Cr, the value of the outstanding notional at T; for simplicity, as 
is commonly done. 

Notice that, differently from what will happen with the tranches 
(see the following section), here the recovery is not considered when 
computing the outstanding notional, in that only the number of 
defaults matters. 

The market quotes the values of So that, for different maturities, 
balance the two legs. If one has a model for the loss and the number 
of defaults, one may impose that the loss and number of defaults 
in the model, when plugged inside the two legs, lead to the same 
risk-neutral expectation (and thus price) 


ee Fol Jo DO, naL | (2.1) 
Bf Eh, DO.) Er] | 


Assuming deterministic default-free interest rates, we can rewrite 


Jp. DO, t)dEof L] 
E & DO, T)(1 — Eo[ Cx ]) 


(2.2) 


So = 
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The assumption of deterministic default-free interest rates will 
turn out to be practical in order to be able to obtain analytical 
or semi-analytical pricing formulas, as is usually done in market 
practice. It can be relaxed to an assumption of independence be- 
tween defaults and default-free interest rates, with the consequence 
that the D(t, T) terms of the deterministic case would become 
P(t, T) = E,[D@, T)], the zero coupon bond prices. 


2.2 CDO TRANCHES 


Synthetic CDOs with maturity T are contracts involving a protec- 
tion buyer, a protection seller and an underlying pool of names. 
They are obtained by putting together a collection of Credit 
Default Swaps (CDSs) with the same maturity on different names, 
1,2,..., M, typically M = 125, each with notional 1/M, and then 
“tranching” the loss of the resulting pool between the points A and 
B,wihO<A<B<1. 


- 1 - 
Le [(LZ; — A)lvcz,<a) + (B — A)liz >B] 
B-A 
An alternative expression that is useful is 
- 1 = S 
La ia [BL?? — AL?] (2.3) 


Once enough names have defaulted and the loss has reached A, 
the count starts. Each time the loss increases, the corresponding 
loss change, rescaled by the tranche thickness B — A, is paid to the 
protection buyer, until maturity arrives or until the total pool loss 
exceeds B, in which case the payments stop. 

The discounted default leg payoff can then be written as 


T 
DerLeca (0) := / D(O, t)dL a? 
0 
Again, one should not be confused by the integral, as the loss ie ki 


changes with discrete jumps. Analogously, the total loss L, and the 
tranche outstanding notional also change with discrete jumps. 
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As usual, in exchange for the protection payments, a premium 
rate sé 2 fixed at time Tọ = O, is paid periodically, say at times 
Ti, D, ..., Ty = T. Part of the premium can be paid at time Tọ = 0 as 
an up-front U 2 The rate is paid on the “survived” average tranche 
notional. If we further assume that payments are made on the notional 
remaining at each payment date 7;, rather than on the average in 
[7;_1, T;], the premium leg can be written as 


PremumLec 4, p(O) — US” + SEP Dva, g(0) 


b 
Dwia,g(0) : = $ & DO, TDO — L7”) 
i=l 
When pricing CDO tranches, one is interested in the premium 
rate S ® that sets to zero the risk-neutral price of the tranche. The 
tranche value is computed taking the (risk-neutral) expectation (in 
t = 0) of the discounted payoff consisting of the difference between 
the default and premium legs above. We obtain 


ia Eo| fo DO, dLE? | E i 
et — (2.4) 
Eo] Dia & DO, TA — Z4”) | 


Assuming deterministic default-free interest rates, we can rewrite 


cp Fa DO, Daal a — ue" 
So = ae 6; D(O, T)(1 = Eo| Li ) (2.5) 


The above expression can be easily recast in terms of the up-front 
premium Us ® for tranches that are quoted in terms of up-front fees. 

The tranches that are quoted on the market refer to standardized 
pools, standardized attachment—detachment points A-B and stan- 
dardized maturities. 

The standardized attachment and detachment are slightly different 
for the CDX.NA.IG and the DJ-iTraxx Europe Main.! 


' The attachment of the CDX tranches are slightly higher, reflecting the average higher perceived 
risk (as measured, for example, by the CDS spread or balance sheet ratios) of the liquid investment 
grade North American names. 
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Table 2.1 Standardized attachment and detachment for the 
DJ-iTraxx Europe Main and CDX.NA.IG tranches 


DJ-iTraxx Europe 


Main CDX.NA.IG 
0-3% 0-3% 
3—6% 3-7% 
6-9% 7-10% 
9-12% 10-15% 
12-22% 15-30% 
22-100% 30-100% 


For the DJ-iTraxx and CDX pools, the equity tranche (A = 0, B = 
3%) is quoted by means of the fair Us B while assuming sé 2 — 500 
bp. The reason for the equity tranche to be quoted as up-front is to 
reduce the counterparty credit risk that the protection seller is facing. 
All other tranches are generally quoted by means of the fair running 
spread sé 2. assuming no up-front fee (Uy 8 — 0), Following the 
recent market turmoil the 3—6% and the 3—7% have also been quoting 
in terms of an up-front amount and a running So 2 — 500 bp given 
the exceptional risk also priced by the market for mezzanine tranches. 


3 


Gaussian Copula Model and 
Implied Correlation 


The Gaussian Copula model is a possible way to model the depen- 
dence of random variables and, in our case, of default times. As the 
default event of a credit reference is a random binary variable, the 
correlation between default events is not an intuitive object to handle. 
We need to focus our attention rather on default times. 

We denote by 1; the default time of name i in a pool of M names. 
The default times of different names need to be connected, and the 
copula formalism allows us to do this in the most general way. 

Indeed, if p;(t) = Q{t; < t } is the default probability of name 
i by time t, we know that the random variable p;(t;) = U; is a 
uniform random variable. Copulas are multivariate distributions on 
uniform random variables. If we call C(u,,...,u,) a multivariate 
uniform distribution, and U,,..., Uy a multivariate uniform with 
distribution C, then a possible multivariate distribution of the default 
times with marginals p; is 


Tı := pī Wy), E TY I= Pu (Um), (3.1) 


where, for simplicity, we are assuming the p’s to be strictly invertible. 
Clearly, since the U;,..., Um variables are connected through a 
multivariate distribution C, we have a dependence structure on the 
default times. 
The Gaussian Copula enters the picture when we assume that 


[U], e...’ Uml] = [®(X)), ONS D(X y)] 


where the X; are standard Gaussian random variables and 
[X1,..., Xy] is a given multivariate Gaussian random variable with 
a given correlation matrix. ® is the cumulative distribution function 
of the one-dimensional standard Gaussian. 
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In the Gaussian Copula model the default times are therefore linked 
via normally distributed latent factors X. 

A particular structure is assumed for the default probabilities p;. 
The default probabilities of single names are supposed to be related 
to hazard rates à. In other terms, p;(t) = 1 — exp(— h ài(s)ds). 

It follows that if t1, ..., ty are the default times of names 1 to M, 
each with a possibly distinct hazard rate A; and cumulated hazard 
rate A;(t) = h Ai(s) ds, then the exponentially distributed default 
times are characterized by 


t = Aj '(—In( — ®(X;))), 
n = Az'(—In(1 — ®(X2))), 


tu = Ay (— Ind — ®(Xm))) 


We also anticipate that the correlation matrix characterizing the 
Gaussian Copula is often taken with all off-diagonal entries equal 
to each other according to a common value p in corporate synthetic 
CDO valuation applications. This is always the case, in particular, 
when defining implied correlations. 

The copula model is used in the market first of all to find a corre- 
lation parameter that is consistent with market quotes for tranches. 
The objects that are quoted by the market are spreads (or sometimes 
up-front prices) of tranches. With the above notation, Se ý (spread) 
or Ue 2 (up front) would be provided by the market and a correlation 
p characterizing a Gaussian Copula with matrix 


1 p p p 
p 1 p p 
p p 1l p (3.2) 
ppp... 1 


would be implied from the market quote. Indeed, Formula (2.4) 
defined the market quotes in terms of expectations of the tranched 
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loss L4-8; in turn, the loss to be tranched at a given time is defined 
in terms of single default times as 


Zyd sE 
L, = 2 ye — les 2 ml T Rd tecnsro 
i= is 
where the X; are multivariate Gaussian variables with correlation 
matrix (3.2) and R; are the recovery rates associated to each name. 
Before getting into more detail on the models, we present two 
general contexts where the copula model can be applied: 


1. Heterogeneous versus Homogeneous Pool. In the homogeneous 
model all single credit references in the CDO pool are assumed 
to have the same probability of default and recovery rates. The 
above p; and R; are all equal to each other. The heterogeneous 
model allows each single credit reference i to have different 
probabilities of default p; and recovery rates R;. 

2. Finite Pool versus Infinite Pool. In the finite pool model we 
take into consideration the actual number of credit references M 
comprising the CDO pool, whereas in the infinite pool model the 
number of credit references is assumed to tend to infinity and 
limit properties are used. 


We will see in particular how combining the assumption of homo- 
geneity with the assumption of infinite number of credit references 
results in a convenient simplification of the tranche NPV calculations. 

The correlation parameter p is therefore affecting the tranche price 
since it is a key statistical parameter contributing to the tranched loss 
distribution, whose expectation will be used in matching the market 
quote. Implied correlation aims at finding the value of consistent 
with a given tranche quote when every other model parameter (the 
pi, the R;) has been fixed. 

We next illustrate the two main types of implied correlation that 
may be obtained from market CDO tranche spreads: 


1. Compound correlation is more consistent at single tranche level. 
However, for some market CDO tranche spreads compound cor- 
relation either cannot be implied or can be implied with multiple 
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solutions, in that there can be more than one correlation value 
reproducing the given tranche spread. 

2. Base correlation is less consistent but more flexible and can be 
implied for a much wider set of CDO tranche market spreads. Fur- 
thermore, base correlation is more easily interpolated and leads 
to the possibility of pricing non-standard detachments. Even so, 
base correlation may lead to negative expected tranche losses, 
thus violating basic no-arbitrage conditions. We illustrate these 
features with numerical examples. 


We will first introduce the general one-factor Gaussian Copula 
model. We will then introduce the finite pool homogeneous one- 
factor Gaussian Copula model. We will show how the loss probabil- 
ities formulas can be computed in this case. The following step will 
introduce heterogeneity, showing that it can be taken into account 
at the expense of efficiency of the calculations, leading to comput- 
ing the loss distribution either by Fourier methods or by recursive 
approaches. We will then go back to the homogeneous assumption 
and show how the infinite pool homogeneous model, popularly re- 
ferred to as LHP (Large Homogeneous Portfolio) makes formulas 
even simpler by means of the law of large numbers. Finally, in this 
work we stay with the homogeneous version, since this is enough to 
highlight some of the key flaws of implied correlation. It has to be 
said that, in practice, due to the necessity of computing hedge rations 
with respect to single names, the heterogeneous version is used, even 
if spread dispersion across names — possibly leading to very different 
pi —complicates matters. 


3.1 ONE-FACTOR GAUSSIAN COPULA MODEL 


A one-factor copula structure is a special case of the Gaussian Copula 
above where 


Xp Sale S +1- iYi, 
with Y;, S standard independent Gaussian variables. S is a systemic 


factor affecting the default times of all names, and Y; is aidiosyncratic 
factor affecting just the ith name. 


Gaussian Copula Model and Implied Correlation 25 


This parameterization would lead to a correlation between Gaus- 
sian factors X; and X; given by ,/p;/P;. However, consistent with 
(3.2) and the homogeneity assumption, all pairwise correlation pa- 
rameters collapse to a single common value, i.e. p; = p for all i. 
We have the same drastic simplification as in (3.2), in that all possi- 
ble 125 x 124/2 = 7750 pairwise correlations across 125 names are 
assumed to be equal to each other. 

Assuming a deterministic and possibly distinct recovery rate R; 
upon default, we are able to simulate the pool loss at any time t 
starting from the simulation of the Gaussian variables X;. 


St 
L, =>) 0 - RDI gen 
i=l M 
M 
1 
= 2 mY — Ri) xico- pO) (3.3) 


i=l 


Thus we are able to simulate the NPV of the premium and default 
leg of any tranche. 


3.1.1 Finite pool homogeneous one-factor Gaussian 
Copula model 


Calculating the NPV of a derivative instrument via simulations can 
be necessary but may lead to intensive numerical effort. Introducing 
the assumption of homogeneity, it turns out that the Gaussian Copula 
model yields a semi-analytical formula to calculate the distribution 
of the pool loss. 

The homogeneity assumption is equivalent to the assumption that 
for each credit reference 


(1) the cumulated hazard rates functions associated with each name 
are all equal to each other: 


Ai(y) =... = AMO) = AQ) 


for all y. Equivalently, this means that the default probability of 
all names in the pool are equal: pı(t) = ... = pm(t) for all t; 
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(2) the recovery rate upon default is the same: R; = R for all i; 
equivalently this means that in order to know the pool loss we 
only need to know how many entities in the pool defaulted and 
the common recovery R. 


Let us assume that we know the realization of the systemic factor 
S. In this case, conditional on S, the default events of the pool of 
credit references are independent. Also in this case, thanks to the one- 
factor assumption and to the homogeneous assumption, we know the 
probability of default of each credit reference in the pool. 

The default by time 7, conditional on the realization of the systemic 
factor S of each single credit reference in the pool, is a Bernoulli 
random variable with the same event probability of default 


-1 = = = 
b (1 — exp(—A(T))) ve"). (3.4) 


Vvl—p 
The number of defaulted entities in the pool by time T conditional 


on the realization of the systemic factor S is the sum of M Bernoulli 
variables and thus is binomially distributed. 


Qi < TIS) = P ( 


! 


= n n —n 
Q{cr=—|s| = JLE < TIS)" -Qr < TIS)” 


M 
ni(M—n 
(3.5) 
where we point out that the Q(t; < T|S) terms are the same for all 
i, since we are dealing with a homogeneous pool assumption. 
We can now integrate the expression in (3.5) to get the uncondi- 
tional probability of n defaults occurring to the standardized pool of 
M credit references before time T. 


afer=H}=[ afër sjeas 6.6) 


—Co 


where ¢(S) is the standard normal density. 


1 62 
-is 


(S) = 


1 
~V 21 
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From (3.6) we immediately obtain the pool loss distribution 
z n > n 
Lr=-R)—} =Qfĉr =<]. 3.7 
Q{Lr=a-R)—} = Q {er = 7 3.7) 


By computing the integral in (3.6) for n = 1, ..., M we obtain the 
unconditional distribution of the pool loss rate that we need in order 
to compute the theoretical tranche spread Ss 2 (or up-front U P) in 
equation (2.4). 

The integral in equation (3.6) does not allow for a closed form 
solution, whereas all other quantities can be analytically calculated. 
Hence the name semi-analytical (analytical up to the calculation of 
the integral) for the finite pool homogeneous one-factor Gaussian 
Copula model formula for the theoretical tranche spread in (2.4). 


3.1.2 Finite pool heterogeneous one-factor Gaussian 
Copula model 


We now remove the homogeneity assumption, and assume that the 
cumulated hazard rates associated with different names can vary: 


Ai(y) Z Aj(y) for i £j 


For simplicity we will retain the homogeneous assumption for re- 
covery rates: R; = R for all i. 

Default of a name i by time T, conditional on the realization of 
the systemic factor S, is a Bernoulli random variable with event 
probability 


bd — exp(—Ai(T))) ve"). (3.8) 


VI=p 

The default counting process, i.e. the sum of these independent 
Bernoulli random variables, is not a binomial random variable, since 
the probabilities are not the same acrossi = 1,..., M. Nevertheless, 
the distribution of the number of defaults (and hence the loss) can 
be obtained via the following algorithm, which basically consists 
of adding one name at a time (one independent Bernoulli random 


Q: < TIS) = o( 
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variable at a time) when building the pool default-counting process, 
all conditional on the realization of the systemic factor S. 

Suppose we know the conditional probability of n defaults occur- 
ring for a subset of K credit references, say names 1, 2,..., K, out 
of the total pool of M credit references, before time T: Q* {Cr = 
n/K|S}, wheren = 1,..., K. QE denotes the probability in a world 
where the pool is constituted only by the first K names. 

Given the conditional event probability of default of name K + 1, 
Q(tx41 < T|S), we can compute the conditional distribution of the 
default counting process for a world with a pool of K + 1 names as 


K+1 JA. _ 0 —_— mk Jr. — 0 
Or = eas] =O fer = RLS] 
x(1 — Q(tk+1 < T|S)),... 


= 2 =] 
QH |e, = zls} = Qk {cr = “= s ecen < TIS) 
kje _ r 
+Q {er = =s} 
x (1— Q(tx41 < TIS), 1<n<K, 
2 K+1 = K 
QH {er = mals = Qk {er = |S} Qtr <T|S). 


At the end of the iterations we would thus get 
Q{Cr = Zis} = Q” {Cr = =s} forall n 
The recursion is initiated by setting 
q! fē = 2 |s =1- Q < TIS) 
[= M ag 1 ’ 
= 1 
1 
Cr=—|s} = Tis 
a {er = FIs] = 26: < r19) 
By computing the integral as in (3.6) for n = 1,..., M, we ob- 


tain again the unconditional distribution of the pool default rate. 
Given the assumption of homogeneous recovery as in (3.7), we 


Gaussian Copula Model and Implied Correlation 29 


obtain the unconditional distribution of the pool loss rate. If re- 
coveries are not homogeneous, the method still works but one needs 
to refine the bucketing of the possible losses for every pool size. This 
method was originally discussed in Andersen et al. (2003). 


3.1.3 Large pool homogeneous one-factor Gaussian 
Copula model 


Let us go back to the homogeneous pool assumption. If we assume 
that the number of credit references comprising the CDO pool tends 
to infinity (or is very large, in practice) then the formulas can be 
simplified further. 

If the portfolio is relatively large, conditional on S the fraction 
of credits of the underlying portfolio defaulting over a specific time 
horizon should be roughly equal to the common individual default 
probability of the underlying credits for that horizon (assumed equal 
across names due to homogeneity). This is a consequence of a prob- 
abilistic result known as the “law of large numbers” (which in turn 
can be derived as a consequence of the Central Limit Theorem). Let 
us see this in more detail. 

Let us introduce the definition of clean pool hazard rate À, taking 
the DJ-iTraxx as a fundamental example: 


DJ-iTraxx DJ-iTraxx 
TAES So = So 


A cican oe = 
1—R Lop 


‘DJ-iTraxx 


where So is the quoted spread of the CDS index for maturity 
T. This formula is based on a single-name CDS analogy where the 
premium leg pays continuously and hazard rates are constant in time 
and deterministic — see, for example, Brigo and Mercurio (2006), 
section 21.3.6. 

If we denote by T the maturity of the contract, then the average 
cumulative default probability of the portfolio can be approximated 
by 


PD(T) = 1 a exp (Aaen T) 
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as in a standard constant hazard rate model. Conditional on the 
systemic factor, S defaults are independent and the probability of 
one default is, as we have seen before, 


o-! (PD(T)) — ves) 
Vvl—p 


Here the Large Pool assumption comes into play: we can assume that 
the number of names going to infinity and check what happens to the 
default fraction in the pool. 

The default fraction, or rate C of the pool at a given time T, is the 
number of defaulted names over the total number of names in the 
pool. We may write, conditional on S: 


PD\(T; S; p) = o( 


1 M 
Cr(S) = M $ leeris 
i=l 


But given that, conditional on S, 1{;,<7)s} are i.i.d. with mean 
PD,(T; S; p), the law of large numbers tells us that 


CHS) —> PD\(T; S; p) as M > œ 


Since, for equal recovery across names, the loss is simply (1 — R) 
times the default fraction C, the loss under infinite pool in our case 
is simply 


L7 (S; p) = (1 — R)PD\(T; S; p). 


If we consider a [0, B] tranche we have that the tranched loss 
(conditional on S) is min (L rS), B). The LHP (see, for example, 
McGinty and Ahluwalia, 2004) computes the expected tranched 
loss by using the above result on the loss under large pool, so that, 
conditional on S, there is no need to take the expectation, due to the 
simplification obtained through the law of large numbers: 


_ z 1 z 
ELES? |S] = LET; S, p) := zo” (L7°(S; p), B) 
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(notice that no expectation is acting really), and, unconditionally, 


BIL} Ko) = IPOT; p) = | E™*™ (7:5, pols) ds 


In particular, substituting all expressions and summing up, we find 


(e (FE) - va 
vl—-p 


7 1 
L°3-(T; p) =f 3min p 


x (1 — R), B | ọ(s)ds (3.9) 


We can value this in closed form, obtaining 


Lop) 


Leo —S™ T 
= O(A)) + z” (-A: 7! ( exp ( E )) ; vP) 


1 — gH T 
a=- oh exp ( 0 )) J1 p 71B /Lo)| 
af. Lop 
and where ®,(-,-,c) is the bivariate standard normal cumulative 
distribution function with correlation c. 
For a generic tranche A, B, we have, through the above formulas, 


L&T; pa, pp) = 


a [BEF on) — A E900) 
where we remember that we are computing everything at maturity 
T. Consistency of the model would impose p4 = pg. We keep the 
possible inconsistency in view of implied correlation. 

Now, let us compute the fair tranche spread. 

Above, we just computed formulas for the expected tranche loss 
for the final maturity T. But since the CDO payout depends on the 
expected tranche loss at different instants, we need to extend it to 
earlier times. McGinty and Ahluwalia (2004) suggest that quarterly 
intervals should be considered, since the payments are made quar- 
terly, and define a “quarterly survival rate” SR in the following way 
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(T is expressed in years, quarterly compounding): 


—4T 
SR? (04, a) 
4 


1 — LAP®(T, pa, pB) = (1 + 


Once SR is backed out for a specific horizon T, we can compute the 
expected notional at any other time T; as 


—4T; 
SR“ (p4, e) 


1 — LP? Tapa) := (1 + i 


Then we can write the tranche default leg price as 


b 


DaLeca,g(0; Pa, OB) = $ DO, Tj) (L^E®(T; pa, pp) 
i=a+1 


—L“®-°(T,_1; P4, Pp) 
and the tranche premium leg for unit spread (or “Dyo1”) as 


b 


Dwa g(0; pa, pa) = >) iDO, TYO — £4? (Tj, pa, pe) 
i=a+1 


In this setup the fair tranche spread sé $ paid in the premium leg, 
which balances the default leg, solves 


DerLec4, g (0; PA, pB) 
Dvo 4, g (0; PA, pB) 


A,B 
So’ = 


This approach has some limits: 


e It can be used only for very large portfolios (typically basket of 100 
names or more). This is the case for CDOs but not for kth-to-default 
swaps. 

e The model cannot produce consistent probabilities for a low num- 
ber of defaults (zero or one). This is due to the fact that the model 
does not consider the actual number of defaults but only the (con- 
tinuous) fraction of the defaulted portfolio. 
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3.2 DOUBLE-t COPULA MODEL 


The Double-t Copula was first introduced — although without much 
detail — by Hull and White (2004). Its main purpose, as we will 
explain in section 3.8 below, is to improve the consistency of implied 
correlation by flattening the base correlation skew, even allowing for 
a flat calibrated skew in some particular dates. Despite the name, the 
reader should pay attention to the fact that the resulting copula is not 
a Student t Copula. 

In the following we will see the details of the calculations only for 
the finite pool homogeneous case, as the extension to the other cases, 
leading to the finite pool heterogeneous and large pool homogeneous 
Double-t Copula, are straightforward. 

Let us consider again the latent factors equation 


Xi = / SH 1 iYi, 


with all Y; independent Student ¢ variables with dı degrees of free- 
dom, and S Student ¢ variables independent of all Y; with d2 degrees 
of freedom. 

Given the homogeneity hypothesis, the cumulated hazard rates 
functions associated with each name are all equal to each other: 


Ai(y) =... = Am(y) = AQ) 


for all y. 
Assume that we know the realization of the systemic factor S. The 
homogeneous conditional probability of default will be: 


F'(A — exp(-A(T))) — VP S 
: 3.10 
T ) (3.10) 


where Fy is the cumulative Student ¢ distribution with d degrees of 
freedom and F, is the cumulative distribution function of the sum 
of two independent Student ¢ distributions with dı and dz degrees of 
freedom respectively (hence a convolution). 

A point that is seldom discussed in Hull and White (2004) is 
that F} is not known in closed form but needs to be computed or 


Ql; = TIS) = Fa ( 


34 Credit Models and the Crisis 


tabulated from the convolution of the two Student ¢ variables. Vrins 
(2009) addresses this task by suggesting either a simulation or a semi- 
analytic approach. We add that one possible way is computing the 
characteristic function of the sum as the product of the characteristic 
functions of the two Student ¢ variables and use inverse Fourier 
methods to obtain the cumulative distribution function F, from the 
characteristic function of the sum. Alternatively, one can use moment 
generating functions and an inverse Laplace transform. 

As in the Gaussian case, conditional on S, the default events of the 
pool of credit references are independent. The number of defaulted 
entities in the pool by time T conditional on the realization of the 
systemic factor S is the sum of M Bernoulli variables, and is thus 
binomially distributed. 


ofer = 2s} 


M! n M-n 
= nim mi oO < TIO A — Q: < T|S)) (3.11) 
where we point out that the Q(t; < T|S) terms are the same for all 
i, since we are dealing with a homogeneous pool assumption. 

We can now integrate the expression in (3.11) to get the uncondi- 
tional probability of n defaults occurring in the standardized pool of 
M credit references before time T. 


z n TS 2 n 
Q{cr ==} = Q{Cr=—|s| Das 619 
where fa, (S) is the Student t density with dz degrees of freedom. 

Thus we see that the only difference in the calculations with respect 
to the Gaussian case lies in the conditional default probability and in 
the integral of the conditional probability of observing n defaults. 

A linear combination of jointly Student ¢ variables is no longer a 
Student ¢, unlike a linear combination of jointly Gaussian variables 
that is still Gaussian. In this respect the above equation for the latent 
factors does not allow us to characterize the one-factor Student t 
Copula. In particular, as we mentioned above, the resulting copula 
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is not a Student t Copula. Nevertheless, it is still a copula whose 
joint distribution we can compute semi-analytically, via a relatively 
simple numerical integration, as we have seen above. 


3.3 COMPOUND CORRELATION AND 
BASE CORRELATION 


Compound correlation is a first paradigm for implying credit default 
dependence from liquid market data. This approach consists of link- 
ing defaults across single names through a Gaussian Copula, where 
all the correlation parameters are collapsed into one. 

Given this correlation parameter, and given the desired declination 
of the one-factor Gaussian Copula, we can compute the loss distri- 
bution on a given set of dates. Thus we can compute the expectations 
contained in (2.4), and with them the fair tranche spread. 

In all calculations below for both the Compound Correlation and 
the Base Correlation, we will make use of the finite pool homoge- 
neous one-factor Gaussian Copula model, avoiding the limit as in the 
LHP model. 

Consider the DJ-iTraxx tranches for example, to clarify the proce- 
dure. For a given maturity in 3,5, 7 and 10 years, consider the market 
quotes 


U®:3% vu 4 500 bp running, 


53-8% er 562% Mi 59-12% ua 51222% Mkt 


To obtain the implied correlation one proceeds as follows. 
First, solve in p3% for the equity tranche 


DirLs09,3%0(0, 03, P3) = 500 bp Dwio,3(0, p3, p3) + UP” 


Then, in moving on one has two choices: retain p3 from the earlier 
tranche calibration and solve 


DerLec3,6(0, 03, p6) = S8% Dvo 6(0, 03, Po) 
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in pe (base correlation), or solve 


) = 52-6% Mkt 


DerLec3, 6(0, 53,6; 3,6 Dvoi3,6(0, 03,6, 03,6) 


in a new £36 (compound correlation). The next step will again be 
similar: for the base correlation solve 


DerLece,9(0, 16, 09) = S°?”™ Dvog, o(0, p6, Po) 


in a new pọ (base correlation) or solve 


§6,9% Mkt 


DerLecg,9(0, 6,9, 6,9) = Dvoi6,9(0, 06,9, 06,9) 


in /6,9(compound correlation). And so on. 

Stop at the 3, 6 tranche for a second: compound correlation is more 
consistent at the level of a single tranche, since we value the whole 
payoff of the tranche premium and default legs with one single copula 
(model) with parameter (3 6. 

The base correlation is inconsistent at the level of a single tranche: 
we value different parts of the same payoff with different models, 
i.e. part of the payoff (involving Lo,3) is valued with a copula in p3, 
while a different part (involving Lo,6) of the same payoff is valued 
with a copula in p6. 

We will focus on implications for base correlation later on, and 
deal now with compound correlation. 

The market data we take as inputs — namely, the reference index 
term structure and 10-year tranche spreads — are detailed in Table 3.1. 


Table 3.1 DJ-iTraxx Europe S5 10-year 
tranches quoted on 3 August 2005 


Tranche Running Up front 


0-3% 500 bp 49% 
3—6% 360 bp 0% 
6-9% 82 bp 0% 
9-12% 46 bp 0% 


12-22% 31 bp 0% 
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Compound Correlation 
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Figure 3.1 Compound correlation smile calibrated to the DJ-iTraxx S5, 10-year 
tranche spreads on 3 August 2005. 


In Figure 3.1 we present the compound correlation smile we imply 
from this set of market data. 

We notice that there is no bar corresponding to the 6-9% tranche: 
from the market spread of the tranche we cannot imply a compound 
correlation. We see in the following how this problem is not atypical, 
in that we quite often face market spreads where we cannot imply 
the compound correlation. 


3.4 EXISTENCE AND NON-MONOTONICITY OF 
MARKET SPREAD AS A FUNCTION OF 
COMPOUND CORRELATION 


We have just seen that on a particular date we cannot imply a com- 
pound correlation from the market spread of the 6-9% tranche. 

To see why this is the case, we investigate this date further, plotting 
in Figure 3.2 the fair market spread as a function of the compound 
correlation: the equity tranche is quoted up front (0.25 means 25%) 
and all other tranches are quoted in number of running basis points 
(123 means 1.23% per annum). The horizontal line is the level of the 
market spread. 

We note that: (1) for certain tranches, from the unique market 
spread we can imply more than one compound correlation, although 
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Figure 3.2 DJ-iTraxx S5 10-year compound correlation invertibility. Tranche 
market spread (horizontal line) versus theoretical tranche spread obtained varying 
the compound correlation between 0 and 1 (dotted line). 


this does not happen in our example of Figure 3.2 where the 
horizontal line crosses the dotted line at most in one point; (2) given 
a market spread we are not always guaranteed that we can imply a 
compound correlation, as we see for example in the 6-9% tranche 
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of Figure 3.2 (there is no intersection between the horizontal line 
and the dotted line). 


3.5 INVERTIBILITY LIMITATIONS OF 
COMPOUND CORRELATION: PRE-CRISIS 


We have seen that on 3 August 2005 we cannot imply a compound 
correlation for the 6-9% tranche. 

In Figure 3.3, taken from Torresetti et al. (2006b), we see how this 
problem is not limited to a sporadic set of dates but is rather affecting 
clusters of dates. 

From March 2005 to November 2006, (the sample on which the 
analysis in Torresetti et al. (2006b) is based), the non-invertibility of 
the compound correlation concerned: 


(1) the 10-year DJ-iTraxx 6-9% tranche (see Figure 7.2 on page 89); 
(2) primarily the 10-year 7-10% tranche and marginally the 10-year 
10-15% CDX tranches (see Figure 7.4 on page 91). 


3.6 BASE CORRELATION 


Here we will illustrate base correlation and will see how compound 
correlation problems are overcome, but at the expense of introducing 
a deeper inconsistency. 


ITRAXX: 0 = compound not invertible CDX: 0 = compound not invertible 
1 1 
04-719 88 1 0 Teen IT 
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Figure 3.3 Compound correlation invertibility indicator (1=invertible, 0=not in- 
vertible) for the DJ-iTraxx and CDX tranches. Source: Torresetti et al. (2006b). 
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Figure 3.4 Base correlation calibrated to the DJ-iTraxx S5 10-year tranche 
spreads on 3 August 2005 (left chart) and expected equity tranche loss (detachment 
set to the standardized tranches detachment) corresponding to the calibrated base 
correlation (right chart). 


In Figure 3.4 we plot the base correlation calibrated to the market 
data in Table 3.1 and the expected equity tranche loss for the various 
detachment points as a function of time. 


EL? ], B = 3%, 6%, 9%, 12%, 22%. 


From these expectations, using equation (2.3), we can compute the 
expected tranche loss E[L A at plotted in Figure 3.5, as a function of 
time. 


From Figure 3.4 we note that the base correlation is a much 
smoother function of detachments than compound correlation. Also, 
to price a non-standard tranche, say a 4—15% tranche, we can interpo- 
late the non-standard attachment 4% and detachment 15%, whereas 
with the compound correlation we do not know exactly what to inter- 
polate (since with compound correlation there is a unique correlation 
associated with each tranche, i.e. correlation is associated with two 
points rather than one). 

Also, as we can see from our examples, the base correlation ap- 
proach is not immune from inconsistencies. In fact in Figure 3.5 we 
note that, in 2005, the expected tranche loss for the 6-9%, 9-12% 
and 12—22% tranches, has already become initially slightly negative. 
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Figure 3.5 Expected tranche loss as a function of time derived from the base 
correlations calibrated to the DJ-iTraxx S5 10-year tranche spreads on 3 August 


2005. 


This inconsistency arises from the different base correlations we use 
in equation (2.3) to compute the two expected tranche loss terms in 
A and B. 
Also, the base correlation still display a considerable skew leading 
to completely different distributions, as can be seen from Figure 3.6. 
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Figure 3.6 Simulated distribution (right tail magnification on the right chart) 
corresponding to a correlation of 15% (dotted line) and 68% (continuous line) 
obtained via a homogeneous finite pool one-factor Gaussian Copula with a 9.3% 
probability of default. 


Let’s consider in fact the 10-year default count distribution corre- 
sponding to a correlation of 15% and 68%. These were the base 
correlations for the detachment 3% and 22% calibrated respectively 
to the DJ-iTraxx S5 10-year tranche spreads on 3 August 2005.! 


3.7 IS BASE CORRELATION A SOLUTION TO 
THE PROBLEMS OF COMPOUND CORRELATION? 


The answer is in the affirmative, and this can be clearly seen for 
example in Figure 3.7 where we plot the fair tranche spread as a 
function of the base correlation on the detachment point for each 
tranche, given the base correlation on the attachment point set equal 
to the calibrated base in the left chart of Figure 3.4. 

This gives us an idea of the range of the tranche spread we can 
calibrate using base correlation. The plots of Figure 3.7 can be com- 
pared with the plots in Figure 3.2, showing the fair tranche spread as 
a function of compound correlation. In Figure 3.7 the thick black line 
is flat at the level of the market spread for the tranche. The two thin 


' The default count distribution is obtained via a homogeneous finite pool one-factor Gaussian 
Copula: the homogeneous 10-year probability of default calibrated from the DJ-iTraxx S5 CDS 
term structure was 9.3%. 
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Figure 3.7 Theoretical tranche spread as a function of the detachment base cor- 
relation keeping fixed the attachment base correlation (dotted line), market tranche 
spread (thick horizontal line) and minimum and maximum tranche spreads invert- 
ible from the compound correlation (thin lines). 


lines are the minimum and maximum spreads we are able to obtain 
by varying compound correlation. 

We note that for each tranche the fair spread is a monotonic function 
of the base correlation on the detachment point, and also that the range 
of market spreads that can be attained by varying base correlation is 
much wider than the corresponding one for compound correlation. 
Consider, for example, the 6-9% tranche: from Figure 3.2 the tranche 
spread that can be inverted in a compound correlation setting lies 
between 93 and 268 bp, whereas from Figure 3.7 the tranche spread 
that can be inverted in a base correlation setting lies in the wider 
range between 0 and 732 bp. 

In Figure 3.7 we did not plot the market spread for all base cor- 
relations between 0 and 1 because, beyond a certain point, the fair 
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Figure 3.8 Expected tranche loss of the DJ-iTraxx S5 10-year tranche for two 
different levels of the detachment base correlation: 38% (left chart) and 48% (right 
chart). 


tranche spread becomes negative. Recall once again that in equa- 
tion (2.3) we use two different correlation parameters for different 
parts of the same payoff: when these two correlations are very differ- 
ent from each other (the detachment correlation is much higher than 
the attachment one) the inconsistency of a negative expected tranche 
loss becomes more evident. 

Consider, for example, the 6-9% tranche. In Figure 3.8 we plot, in 
the abscissas, the year fraction of the tranche payment dates and, in 
the ordinates, the expected tranche (6-9%) loss. For both graphs in 
Figure 3.8 the tranche attachment correlation is the calibrated base 
on the 6% detachment. The tranche detachment (9%) correlation is 
set to the calibrated base on the left-hand graph (38.07%) and to an 
arbitrarily high level (48%) on the right-hand graph. 

We can see a markedly negative profile for the expected tranche 
loss, which clearly violates no arbitrage constraints. Indeed, the loss 
must be non-negative and non-decreasing in each path through time, 
and tranching does not alter that. More simply, the tranche loss at 
a given point in time is a non-negative random variable, and its 
expected value also needs to be non-negative. This does not happen 
in our example as the basic non-negativity constraint is violated. This 
is a strong disadvantage of the base correlation paradigm at a very 
basic level. 
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3.8 CAN THE DOUBLE-t COPULA FLATTEN THE 
GAUSSIAN BASE CORRELATION SKEW? 


We have already mentioned at the beginning of section 3.2 that 
the Double-t Copula was introduced in order to try to flatten the 
pronounced base correlation skew displayed by the Gaussian Copula. 
Finding a flat implied base correlation curve along the detachment 
dimension means that the tranches are all reproduced with the same 
parameter. This would give consistency, since now one would have a 
single Double-t Copula that is consistent with all the tranches at the 
same time for a given maturity. 

Previous analysis of Hull and White (2004) shows that these results 
occurred with the Double-t Copula in a particular date for dı = d? = 
4. They present an example based on the DJ-iTraxx on 4 August 
2004, where the implied correlation across tranches ranges between 
0.24 and 0.266 using a Double-t specification, whereas it ranges 
between 0.204 and 0.448 with the base correlation obtained from the 
Gaussian Copula. 

Even if the Double-t Copula is seen to be better in this case, it is 
hardly to be expected that during the crisis — especially with the very 
high values taken by super-senior tranches — the Double-t Copula 
could fit the data with a flat base correlation curve across detach- 
ments. It might be that pairing the Double-t Copula with random 
recovery can help, but the framework becomes numerically intensive 
and possibly unstable due to the need to numerically evaluate F}. 


3.9 SUMMARY ON IMPLIED CORRELATION 


Calibrating base correlations by decomposing the tranched loss, as 
shown in equation (2.3), solves the two main issues concerning com- 
pound correlation calibration. 

Indeed, as illustrated in Figures 3.4 and 3.7, by calibrating directly 
base correlations: 


(1) we get a monotonic mapping of the base correlation parameter 
into tranche spreads; 
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(2) we can invert a wider range of tranche spreads into a base corre- 
lation parameter; 

(3) we can price bespoke detachments by interpolating the base cor- 
relation across detachments. 


Even so, base correlation needs to be handled with care, since it may 
lead to negative losses, thus violating basic no-arbitrage constraints. 

Also, as with compound correlation, 125 x 124/2 pairwise depen- 
dence parameters are compressed into one. Assuming a flat correla- 
tion across the matrix, characterizing the copula, is a gross simplifi- 
cation that does not reflect reality at all. 

Alternative copula specifications are possible. Indeed, Hull and 
White (2004) show that on a particular date the Double-t Copula can 
consistently reproduce tranche spreads without skew in the correla- 
tion parameter. Our independent tests show that the skew resurfaces 
at later dates (from May 2005 onwards). But the key problem of 
the copula — being a static concept — cannot be solved by changing 
distribution. 

A more model independent approach to tranche interpolation and 
pricing consists of implying expected tranche losses without assum- 
ing any model (see, for example, Walker, 2006; Torresetti et al., 
2006a; and Cont and Minca, 2008). 

There are hence two major problems with the use of implied cor- 
relation as a parameter in a Gaussian Copula model, as we described 
above, even if one is willing to accept the flattening of 7750 param- 
eters into one. First, the model parameter changes every time we 
change tranche, implying very different and inconsistent loss distri- 
butions on the same pool. Second, there is a total lack of dynamics 
in the notion of a copula. We have addressed the first problem here, 
showing that the financial literature has been doing that well before 
the current crisis started in 2007. We will address the second problem 
below. 


4 


Consistency across Capital 


Structure: Implied Copula ___| 


In the Implied Copula approach, a factor copula structure is as- 
sumed, similar to the one-factor Gaussian Copula approach seen 
earlier. However, this time we do not model the copula explicitly, 
but we model default probabilities conditional on the systemic fac- 
tor S of the copula: the copula will then be “hidden” inside these 
conditional probabilities that will be calibrated to the market. Hence 
the name “Implied” Copula. In illustrating the Implied Copula we 
will also assume a large pool homogeneous model in that the default 
probabilities of single names will all be taken equal to each other. 

Let us consider, for simplicity, survival probabilities that are as- 
sociated to a constant-in-time hazard rate. We know that if we have 
a constant-in-time (possibly random) hazard rate à for name i, then 
the survival probability is 


Q(t > t) = Efexp(—Ar)] 


The Implied Copula approach postulates the following “scenario” 
distribution for the hazard rate A conditional on the systemic factor S: 


Conditional Systemic Scenario 
hazard rate scenario probability 
Ài S= sj pı 
A|S ~ n 
À2 S=s p2 
Àn S = Sp Ph 
This way the default probability for each single name i = 1,..., M 


is, conditional on the systemic factor S, 


Q(t; < t|S =s5;) = 1 —exp(—Ajt). 
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Compared with the Gaussian factor copula case: 


-I(l — exp(—A(T))) — ves) 
Vl—p 


Unconditionally, the Implied Copula yields the default 
probabilities 


Q(z; < ris=s) = 


h h 
Qa < t) = $ piQ < tIS = sj) = Ñ p; (1 — exp(—Ajt)) 


j=l j=l 


Conditional on S, all default times are independent, have the same 
hazard rate and their hazard rates are given by the above scenarios. 

By resorting to the infinite pool approximation, as we illustrated for 
the Gaussian Copula in the LHP version, we have that the conditional 
default indicators 1{7,<7r|s=s,) are i.i.d. for all 7, so that their sample 
average as their number tends to infinity tends to the single common 
true mean: 


Cy (sj): a> Lt,<7|s=s;) > Ele, <r|s=s,) 
Su ihe 


when M tends to oo. 

Again, by this method we avoid taking expectations, except the 
final one with respect to S, since, conditional on S = sj, all random- 
ness has been ruled out by the law of large numbers and both the 
default rate and the loss are completely determined. 


Systemic Scenario 
Conditional default rate scenario probability 
= CM (si) = 1 — exp(—A,T) S=s | 
M A ~T 
Cr) ) EMs) = 1—exp(-aaT) S=s p SD 


CM (sn) = 1—exp(—AnT) S=Sn Ph 
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We will assume, in line with the market convention, that the protec- 
tion payment will be calculated on the average outstanding notional 
between any two protection premium payment dates. 


Prenuunilc’y g T= Cat + So’? Dvo} g 
. b . 
Dwi, g 2 =). 5 DO, T; — L(0.5(T; + Ti-1))4 8) 
i=l 
j min(B — A, max((1 — R;)(1— exp(—A j;t)) — A, 0)) 
L(t), pi = 
i B-A 


where U Band sé 2 are the market mid up-front and running spread 
for the tranche A, B. 

We will discretize the loss increments, entering the calculation of 
the discounted default leg payoff, on the same set of dates of the 
discounted premium leg payoff calculation: the protection premium 
payment dates. We will also assume that, on average, the loss incre- 
ment arrives at the middle of each time interval [7;_,, 7;]. 


b 
DerLeok g : = (1 — Rj) X D (0, 0.517; + Ti-1)) (LUT, g 
i=l 
-L(T;-1)4 g) 
where R is the deterministic recovery rate typically set to 40% for 
Senior Unsecured CDSs: the same type of CDSs, comprising the 
DJ-iTraxx and CDX standardized pools. 
We will call PremumLeca,g and DerLecg g the column vector that 


stacks respectively all discounted premium leg and default leg values 
conditional on the systemic factor s states: 


PremumLec A,B >= [Premec B œ> PremumLect Pale 
DeL = [DeLec! Dileh 31" 
EFLEGA,B >= ELEGA Bootes EFLEG 4 Bl - 
The receiver tranche price conditional on S = s; is simply 


Now, B = PremumLesh B E. DerLec’y B 
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Then we integrate against S, simply summing over all possible haz- 
ard rate scenarios multiplying by the scenario probability, to obtain 
the tranche unconditional price: 


h 


Nev 4B = ) PiN g 
j=l 


In matrix notation the receiver tranche value can be rewritten as 
Nev4,B = PT (Premumleo 4, g = DerLec 4B) 


where P = [pı ... pp]? is the column vector with the systemic 
factor probability distribution. 


4.1 CALIBRATION OF IMPLIED COPULA 


When calibrating all the 5-year DJ-iTraxx tranches we end up mini- 
mizing a constrained sum of squares. If we call Nev the matrix with 
all tranches (columns) discounted, the payoff for all possible states 
(rows) is: 
1 1 
Nevo, 3% -+ Nev22%,100% 
Nev = 
h h 
Nevog,3% +++ Nev22%, 100% 
Then, calibration becomes simply 
argmin P7 Ne Nev! P (4.2) 
P=(p1;-- Ph)" 


subject to 


h 
paf p20 i=l1,...,h 
i=l 


Note that for the base correlation calibration the market tranches 
had to have consecutive adjacent attachment and detachment points 
spanning the whole capital structure; if that was not the case, some 
base correlation interpolation assumption had to be introduced. With 
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the Implied Copula this is not necessary in that the tranches to be 
calibrated need not span the whole capital structure or be adjacent. 

In the Implied Copula framework, Hull and White (2006) cali- 
brate the scenario probabilities p while pre-assigning the hazard rate 
scenarios À; exogenously. One may also calibrate the index itself 
besides the tranches, by adding the squared difference between the 
model implied index and the market index in the target function. 

One may naively think that since we have seven quotes (six 
tranches plus the index) then seven parameters p could be enough to 
fit the market quotes, getting a low squared error out of the minimiza- 
tion. So one could set the number of scenarios to h = 7. Actually, 
this does not work. One may have to go up to h = 30 in order to be 
able to fit market data with a good precision. 

We tried the model and checked that, indeed, even with 27 points 
the error remains large, no matter how smartly one preselects the 
Aj, +++, A27. However, even with h = 30 one more problem surfaces. 
If we do not impose some further constraint or smoothing on the 
minimization problem, we are finding multiple solutions each time 
we change the initial guess. We will address this below through a reg- 
ularization technique, but first we discuss some recovery modelling 
issues. 

In this respect, we mention a relationship between default prob- 
abilities and recovery rates that may be necessary to fit the market 
correlation skew in periods of turmoil (e.g July 2005). 

Following the results of an empirical study by Hamilton et al. 
(2005), Hull and White suggest a recovery change in each scenario 
by linking it to the conditional probability of default in that scenario: 


R; = 0.52 — 6.9(1 — exp(—Aj;5y)). 


This functional form will assign a recovery of 0% to all default 
rates above 8%. Alternatively one could fit a logarithmic relationship 
between recovery rates and default rates: 


R; =a log(1 — exp(—Aj5y)). 


52 Credit Models and the Crisis 


Remark 4.1.1 (Random recovery as a function of the systemic 
factor). The above formula makes the recovery a function of the 
intensity conditional on the systemic factor. This way recovery be- 
comes a function of the systemic factor. An approach that makes 
recovery a function of the systemic factor in the inconsistent base 
correlation framework is given in Amraoui and Hitier (2008). 


The dots in Figures 4.1 are the data points we have used to fit the 
parameter a: issuer-weighted all-rated bond default rate and issuer- 
weighted recovery rate senior unsecured from 1982 to 2004 taken 
from Moody’s. From 1982 to 2004 the issuer-weighted all-rated bond 
default rate has always been below 3.82%. Therefore, in Figure 4.1 
we cannot fully appreciate the different implication of the two func- 
tional forms when applied to the Implied Copula. 

In Figure 4.2 we present the two different recovery functions on the 
entire range of possible default rates: from 0% to 100%. We note that 
for a default rate above 8% the recovery obtained via the functional 
form suggested by Hamilton et al. (2005) is 0%. 

Once we select our favorite functional relationship between re- 
covery rates and default rate, the formulas above remain unchanged, 


0% 1% 2% 3% 4% 


+ RR 
Logarithmic Fit 
- - - Hamilton et al (2005) 


Figure 4.1  Issuer-weighted all-rated bond default rate versus issuer-weighted 
recovery rate senior unsecured. Data from 1982 to 2004 (Source: Moody’s). 
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0% E TER "= T 
0% 20% 40% 60% 80% 100% 


Logarithmic Fit 
- - - - Hamilton et al (2005) 


Figure 4.2 Different functional forms between default rates and recovery rates: 
logarithmic (solid line) and floored linear (dotted line). 


except for the discounted default leg payoff conditional on S = sj: 


DeLecoh p: = X D0,0.5(T; + T-1)) 
i=l 
x (exp(—à;Ti—z1) — exp(—àÀ;TiDA — R;) 


Going back to regularization, Hull and White (2006) suggest 
adding to the target function a quantity that penalizes changes in 
convexity in the patterns of the scenario probabilities plotted against 
the default probabilities associated to each scenario: 


(Pj + pj-1 — 2p; 
0.5 (e —Aj-1T, _ eit) 


argmin PT NwNw! P+c ae (4.3) 


P=(p1,---, Pn)? 


subject to: 


a1: POU i=1,...,h 


i=l 
The new term’s numerators are second-order differentials of the 
p terms, and they penalize departures from linearity. Therefore, 
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when we plot the graph {(i, pj), i = 1,..., h} all kinds of humps 
are penalized. The constant c controls the degree of smoothing: too 
large a c privileges the smoothing over the fit, and we cannot recover 
exactly the market quotes; while too small a c fails to guarantee 
substantial uniqueness in the calibration solution. We found good 
results with c = 1. 


4.2 TWO-STAGE REGULARIZATION 


Torresetti et al. (2006c) propose a different approach with respect to 
Hull and White (2006), where there is no need to select a regulariza- 
tion coefficient. 

The main differences between the two approaches are summarized 
as follows: 


e Assign 125 possible states to the system factor S = s;. The hazard 
rates associated with each state are this time such that the pool 
default rate at maturity is equal to (j — 1)/125: 


A; such that (j — 1)/125 = 1 — exp(—Tpà j) 
where j = 1,..., 125. (4.4) 


e In our version perform a two-stage optimization that will assure 
that all tranches are priced within the bid-ask spread without the 
need to choose a regularization coefficient c. 


The numerical solution P of the first stage optimization, i.e. the 
solution to the problem set out in (4.2), can be expressed as a deviation 
of each tranche theoretical spread from the market mid quote. If 
we divide this difference by half the bid-ask spread we obtain the 
standardized mispricing: 


A,B 
So noe Neva, p/Dwi4,B 


GSS ASE = gaa bid 


MiseSmz4, B = 
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where Dw, g is the expected value of the basis point value of the 
tranche 
125 
Dvoi4,B = X pjDvork g 
j=l 
Values of A and B are given by the canonical attachment and 
detachment points of the tranches for the pool under consideration. 
If the standardized mispricing is between —1 and +1, then the 
vector of systemic state probabilities, P, results in a calibration of 
the theoretical tranche spread Nev4,g /Dvo4, that lies within the bid- 
ask, spread. 
The two-stage optimization consists in the following: 


e The first stage consists of the optimization set out in equation (4.2). 
Assuming that the solution results in a theoretical tranche spread 
that is within the bid-ask spread, then we pass to the second stage. 

e This stage consists in seeking the most regular distribution while 
still pricing all tranches within the bid-ask spread. The numerical 
optimization starts from the solution found at the end of the first- 
stage optimization. 


The solution to the first-stage optimization gives a numerical start- 
ing point to the second-stage optimization problem: 


h— 
Hatan 
0.5/125 


argmin 
P=(p1,.-Pn)? j=2 
subject to: 


h 
J pai p20 Fath 


i=1 


MiS? < 1 


Note that in the definition of MiserSwz4’? both Neva.p and Dvoi4 g 
depend on the vector of systemic state probabilities with respect to 
which we are optimizing. 
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Table 4.1 CDX.NA.IG 5-year tranches quotes on 6 June 2006 


Mid Bid Ask 
0-3% 32.40% 31.50% 33.60% 
3-7% 106.5 bp 100 bp 113 bp 
7-10% 23.5 bp 22 bp 25 bp 

10-15% 10 bp 9.3 bp 10.7 bp 

15-30% 5.5 bp 5.1 bp 5.9 bp 


We now put both approaches to work, the one of Hull and White 
(2006) and that of Torresetti et al. (2006c), on a specific date and we 
highlight their differences. 

To make the two approaches comparable we have assigned the 
same 125 possible states to the system factor S = s; in both ap- 
proaches to the Implied Copula. The hazard rate associated with 
each state are such that the pool default rate at maturity is equal to 
(j — 1)/125, as in (4.4). 

We then solve the optimization problem set out in (4.3). We con- 
sider the market quotes in Table 4.1 for the CDX.NA.IG 5-year 
tranches on 6 June 2006. 

In Figure 4.3 we show how varying the regularization coefficient 
c in (4.3) affects the pricing accuracy of the tranches. In the abscissa 
of the charts log), c is shown, so that the regularization coefficient c 
in the charts varies between 1075 and 10°. The two flat black lines 
represent the bid-ask spread for each tranche. In other words, it is 
the maximum positive and negative Nev4,g such that the standardized 
mispricing MiseSmz^? is still smaller than 1 in absolute value. 

In Figure 4.3, on opposite page, we note that if we were to price all 
tranches within the bid-ask spread we would need to choose c = 10!. 
In fact we see that for c > 10! the mispricing of the 7-10% tranche 
falls outside the two flat black lines. 

In Figure 4.4 instead we look at the optimal distributions solv- 
ing the problem set out in (4.3) when the regularization coefficient 
varies between 1075 and 10°. We note that the distribution changes 
considerably in the various cases. 
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Figure 4.3 Tranche mispricing as a function of the regularization parameter c. 
In the abscissa, log), c is shown, i.e. 5 reads c = 10°. 


While in the optimization outlined in (4.3) we do not control di- 
rectly for the tranche mispricing, this could be controlled iteratively: 
find the largest regularization coefficient c such that all tranches are 
priced within the bid-ask spread. 
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Figure 4.4 Optimal pool default counting MCs), distribution obtained varying 
the regularization coefficient c. 


In Figure 4.5, on following page, we plot the optimal distribution 
resulting from the two-stage optimization as outlined in Torresetti 
et al. (2006c) . We see that it is very much in line with the optimal 
distribution that would be obtained in (4.3) setting c = 10!, as can 
be seen by comparison with Figure 4.4. 

This is also true for the tail of the distribution, as can be noticed 
comparing Figure 4.5 and Figure 4.6, on following page. 
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Figure 4.5 Optimal pool default counting MCs,, distribution resulting from the 
two-stage optimization. The right-hand chart magnifies the distribution tail. 
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Figure 4.6 Magnification of the tail of the optimal pool default counting MC 5yr 
distribution obtained by varying the regularization coefficient c. 


4.3 SUMMARY OF CONSIDERATIONS AROUND 
IMPLIED COPULA 


This is a powerful approach with some strong limits. On the positive 
side, we have an approach that is capable of calibrating more CDO 
tranches consistently across the attachment dimension, improving 
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on the lack of consistency of the implied correlation. A drawback 
of this approach, like all other copula approaches, is that the core 
dependence across defaults is modelled as static, and there are no 
dependence dynamics. Pricing a tranche option would require, in 
principle, a dynamical model. 

Also, there is no differentiation for single-name probabilities, 
which are the same (although scenario based) across names. 

The Implied Copula as presented above can solve the problem 
of consistency with tranche quotes for a single maturity: there is a 
unique implied loss distribution that is consistent with all the avail- 
able tranche quotes for a single maturity. Let us assume that one is 
momentarily willing to leave aside the problem of single-name mod- 
elling and hedges, which is implicitly neglected by the homogeneous 
pool assumptions. 

This may be momentarily acceptable, based on the remark that the 
no-arbitrage expressions for the quoted spread of credit indices, given 
by equation (2.1), and of CDO tranches, given by equation (2.4), 
show that the only information we can infer from these market quotes 
are expected aggregated quantities, while we lack direct information 
about dependencies across single names. In particular, credit indices 
depend both on the expected portfolio cumulated loss and on the 
expected number of defaults, while CDO tranches depend only on 
the expected tranched portfolio cumulated loss. 

We can momentarily leave single names aside and focus on the 
fundamental issue of having consistency of the implied loss distribu- 
tions at a single time. However, when CDO tranche quotes for more 
than one maturity are available, we cannot account consistently for 
them with an Implied Copula model. We need a arbitrage-free dy- 
namic loss model, implying future loss distributions for the loss at 
different times that are consistent with the tranche quotes at each 
maturity. 

An explicit model implying a dynamic for dependence across de- 
faults, absent in the copula case, is the above mentioned GPL model 
given in Brigo et al. (2006a, 2006b), and to that we turn now, leaving 
aside the single-name consistent but combinatorially complex GPCL 
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model in Brigo et al. (2007). Finally, we note that the loss distribution 
of the pool under the risk-neutral measure, which has been discussed 
in this book, is different from the actual loss distribution in the ob- 
jective measure. This has been pointed out, for example, in Torresetti 
et al. (2006c), who find the well-known presence of a risk premium 
explaining this difference. See also the discussion in Torresetti and 
Pallavicini (2007) concerning CPDOs, where the difference in the 
modes of the tail loss distribution is also visible. 


5 


Consistency across Capital 
Structure and Maturities: 


Expected Tranche Loss 


In this chapter we approach the issue of consistently representing the 
information embedded in tranches across attachments and maturities 
in a model-free way. 

Let us again consider equation (2.5). 


fË DO.nd Eo| LA? | =i" 


A,B 


° 5, DO, 7) (1 = Eo| 24” |) 


The numerator and denominator of the above formula depend lin- 
early on the expected tranche loss Eo[L - B] at different times (ETL).! 


If a term structure of tranche upfronts, and spreads S a for differ- 
ent maturities are given, then it is possible to strip back the expec- 
tations Rol Ly" ] in a model independent way. We will address this 
issue in this chapter. 

We first notice the nested structure of the expected tranche losses 
in time. Indeed, the NPV calculations of two tranches with the same 
attachment and detachment points (A, B), but different maturities T; 
and T>, both depend on the expected tranche loss Eo[L al for each 
T; < min{7|, T2}. For example, in the NPV calculations of a 5-year 
maturity 0-3% tranche and a 10-year maturity 0-3% tranche, the 
quantity Rol Lz” | should be the same for the calculation of the 5- 
and 10-year tranches for all T; < 5 years. 


' Let us stress again that the actual outstanding notional in each period would be a daily average 
over [T;—1, T;], but for simplicity we replace it with 1 — Lr s the value of the tranche outstanding 


notional at T;. 
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= 3yr -e 5yr -—10 yr 


Figure 5.1 Expected equity (0-3%) tranche loss resulting from the base corre- 
lation, One-factor homogeneous finite pool Gaussian Copula, calibration to the 3-, 
5- and 10-year 0-3% CDX tranches on 26 April 2006. 


In the copula framework unfortunately we are not guaranteed that 
this is the case, regardless of the particular copula model considered, 
Gaussian versus Double-t, or Homogeneous versus Heterogeneous. 
This has to do with the copula being a static object. As an illustration, 
in Figure 5.1 we plot the expected 0-3% tranche loss resulting from 
the base correlation, one-factor homogeneous finite pool Gaussian 
Copula, calibrated to the 5- and 10-year 0-3% tranches on 30 January 
2009. 

We note that the expected 0-3% tranche loss resulting from the 
calibration to the 3-, 5- and 10-year 0-3% tranches do not overlap. 
Thus when valuing the same expected tranche loss Eo Ly] for 
T; < 3 years we are using three different numbers depending on the 
tranche maturity, even though the pool of underlying credit references 
is the same. 

In the following we will see the details of a non-parametric model- 
free estimation of the basic quantities entering linearly the tranche 
NPV calculation, namely the expected tranche losses Rol Ly” ; 
which is consistent through maturities and attachments. 
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5.1 INDEX AND TRANCHE NPV AS A FUNCTION 
OF ETL 


Let us assume that we are given the market spreads of a set of tranches 
spanning the entire capital structure. In particular, let us assume that 
we have the market spreads of a set of k tranches with attachments 
A; and detachments B; with j = 1,...,k, where A; = 0, B = 1 
and Aj; = B; fori = 1,...,k —1. 

Let us consider the standardized DJ-iTraxx tranches with de- 
tachments (3%, 6%, 9%, 12%, 22%): all except the super-senior 
22-100%.* 

Our goal is to price the 5-year tranches and the index: a total of 
six market quotes. To achieve this goal we will be looking for the six 
unknown 5-year ETLs that set the NPV of the instruments as close 
as possible to zero. 

To compute the NPV of the tranches and index, we also need the 
ETL on all payment dates with maturity 7; shorter than 5 years: 
these will be obtained interpolating for each tranche between the 
initial time 0, where by definition the ETL will be zero,’ and the 
six 5-year unknown ETLs that are being sought: 3%, 6%, 9%, 12%, 
22% and 100%. Therefore, once we set a possible value for the 
unknown 5-year ETLs in the iteration implementing the optimization, 
through interpolation we also set a value for all the ETL paths up 
to 5 years that is consistent with the selected 5-year value being 
iterated. 

Once the 5-year nodal ETLs matching the 5-year data are found, 
to price the 10-year tranches and index we will also need the ETLs 
between 5 and 10 years. As before, for each tranche this will be 


2 At the time when Walker (2006) and Torresetti et al. (2006a) ran their analysis, the super-senior 
tranche was just beginning to have some liquidity so that it was not considered in their historical 
calibration analysis of the expected equity tranche loss surface. In fact, before the end of 2005 the 
liquidity of the super-senior tranche was quite limited. 

3 In case the pool of underlying credit references already suffered losses, then a new version 
of the index is rolled out. The existing positions in both the index and the tranches are usually 
promptly rolled in the new version. Tranches of the newly rolled version see the standardized 
attachment and detachment adjusted to reflect the loss in subordination and the reduction in the 
pool outstanding balance. 
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obtained by interpolation between the expected tranche loss at 5 
years and the six unknown ETLs at 10 years to be found. 

We call f(t, h, k) the ETL at time t of the tranche with attachment 
h and detachment k. To simplify the notation we will often identify 
the seniority of the tranche in the capital structure of the CDO only 
through the detachment point k, writing f(t, h, k) = f(t, k) when h 
is clear from the context. 


f(t, A, B) = Eo[ £47 J. 


By assuming interest rates and default times to be independent, 
and discretizing the integrals entering the NPV of the premium and 
default leg in equation (2.4) on the tranche payment dates, we obtain 


s42 _ Lier DO, S(T, A, B) ~ f(T, A, B) 
ú X, & DO, TOA — f(T, A, B)) 
Tj-1 +T; 


ti = a, (5.1) 


Then, since the tranches are adjacent, the expected portfolio loss is the 
summation of the ETLs multiplied by the tranche depth (detachment 
minus attachment), namely: 


k 
f(t, 0, 100%) = SELL" ](B;- A) = EIL] (65.2 
i=1 


Once given the expected portfolio loss E[L;] and the recovery rate 
R, we can compute the expected portfolio default rate as E[L,]/( — 
R) = f(t, 0, 100%)/(1 — R). Hence, the index spread becomes: 


go.to0% _ Liter DO, HFT, 0, 100%) — f(Ti-1, 0, 100%)) 
= X, DO, TAA — f(T, 0, 100%)/(1 — R)) 
(5 


We will seek to calibrate tranches and indices for the 3-, 5-, 7- and 
10-year maturities in terms of the unknown f(t, A, B), which are 
listed on next page: 
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0-3% 3—6% 6-9% 9-12% 12-22% 22-100% 
t=0 0 0 0 0 0 0 
LaS fyr, 3%) fGyr,6%)  ... f Gyr, 100%) 
re 5 fG6yr,3%) fGyr,6%)  ... f Gyr, 100%) 
t=7 S(7yr, 3%) f(Tyr, 6%) ... f(Tyr, 100%) 
t = 10 fayr, 3%) f(10yr,6%) ... f (0yr, 100%) 


The set of 4 x 10 x 6 = 240 f’s (one for each quarterly payment 
date and for each tranche) created by interpolating the 4 x 6 = 24 
basic nodal f(t, k)’s (one for each maturity and for each detachment) 
will be used to set the NPV of the instruments as close to zero as 
possible while maintaining the following constraints 


argmin > > MiseS mz? 4. p (5.4) 


( f Byr.3%),...5 f U0yr, 100%)} F AB 
subject to: 


0< ft,k)<1l 
f(ti,k) = f(ti-1,k) (5.5) 
f(t, kj-1) 2. f(t, kj) 


where 


A,B A,B,mid 
ST, 7 ST, 


A,B,ask  ¢A,B,bid 
(S7, aa — S7, *)/2 


(5.6) 


MisprStvzr, A,B = 


and the double summation in the objective function* is taken with 
respect to all four maturities (3, 5, 7 and 10-year) and all available 


4 Our objective function to be minimized will be the sum of the squared standardized mispricings 
in equation (5.6). The index standardized mispricing will be multiplied by 2 in order to facilitate a 
much closer repricing of the index. If using the sum of the squared mispricings will yield a solution 
for which some of the instruments are priced outside the bid-ask spread (|MisprStpz7,4,2| > 1) 
then we will try to minimize the sum of even powers of the standardized mispricing in (5.6) for 
exponents 4, 6 and 8. In the case of persistence of instruments priced outside the bid-ask spread, 
we take the solution for which the maximum absolute standardized mispricing is smallest. 
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instruments, tranches plus the index (the latter being identified by 
the attachment of 0% and detachment of 100%). 


(A, B) € {(0, 3%), (3%, 6%), (6%, 9%), (9%, 12%), (12%, 22%), 
(0%, 100%)} 


Note that the ETLs f(t, k) are entering the equation of the standard- 
ized mispricing (5.6) of the tranches and indices a B via equations 
(5.1) and (5.3) respectively. 

In Figure 5.2, shown on following page, we plot the ETL f(t, k) 
for a series of dates. We note that the surface is not particularly 
smooth. To price a tranche on the same pool of credit references with 
non-standard attachment and detachment A and B we would need 
to obtain the expected equity tranche losses g(T, B) = Bf(T, 0, B) 
and (T, A) = Af(T, 0, A). 

The optimization problem set out in equation (5.4) in terms of ETL, 
f(t, k) can be formulated in terms of the much smoother Expected 
Equity Tranche Loss (EETL) surface, g(7, B) = Bf(T, 0, B). 


g(t, B) = E[L®?] = Efmin(L,, B)] (5.7) 


Thus, now the problem will be to look for the 12 unknown EETLs 
setting to 0 the NPV of the 12 market instruments: five tranches and 
an index for both the 5- and 10-year maturities. In this approach the 
tranches are adjacent, so that from the EETLs we can always retrieve 
the ETL needed to price the tranche, given that 


f(t, A, B) = B E[min(L,, B)] — A E[min(L,, A)] 
B-A 


E B-A 


(5.8) 


We note that g(t, 1) = E[L,]. While f(T, A, B) represents the ex- 
pected tranche loss and as such is bounded between 0 and 1, g(7, B) 
is bounded between 0 and the expected pool loss E[L;]. 
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Figure 5.2 ETL calibrated to the Dj-iTraxx 5-year and 10-year tranches plus 


index. 


The constraints (5.5) of the calibration in equation (5.4) will be 
restated as follows in terms of (T, k): 


O<gt,k) <1 


g(t;, k) = g(ti-1,k) 
g(t, kj-1) < g(t, kj) 


g(t, kj+1) < g(t, kj) + (kj+1 


(5.9) 
g(t, kj) > a(t, kj-1) 


kj) 
f kj —kj-1 
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The first and second constraints specify that, since the pool loss 
paths are bounded within [0, 1] and are non-decreasing, the EETL 
must do the same. The third constraint requires that, by increasing the 
cap B in equation (5.7), expectations with respect to the detachment 
must be non-decreasing. The fourth constraint is simply equivalent to 
requiring that the ETL f(T, h, k) is decreasing with respect to k, see 
equation (5.8), so that the constraint that f(t, k) must be decreasing 
in k is equivalent to requiring that g(t, k) must be convex in k. In 
turn, this is equivalent to requiring its derivative with respect to k to 
be decreasing (or its second derivative to be non-positive). 


F(t, Kj-2, &j-1) — f(t, kj-1, kj) 
ee kj- f(t, 0, kj-1) = kj f(t, 0, kj-2) 


he i kj-2 
kj f(t,0,kj) —kj-1 f(t, 0, kj-1) T 
k; — ki 
an g(t, kj-1) — g(t, kj-2) x g(t, kj) — g(t, kj-1) 
ae oe ne 


so that 
g(t, kj) — g(t, kj-1) z 
kj = kj-1 7 
Thus requiring a decreasing ETL f(t, k) with respect to k, i.e. 
as we move on the capital structure towards larger detachments, is 
equivalent to requiring a decreasing derivative of the EETL g(t, k) 
with respect to k. This in turn is equivalent to requiring that the EETL 
be convex in k. This can be seen more rigorously through a Breeden 
and Litzenberger (1978) argument, by differentiating twice the EETL 
with respect to the detachment. This leads to a quantity proportional 
to the opposite of the risk-neutral density of the loss distribution 
computed at the detachment, which is clearly non-positive. 
Comparing Figure 5.2 to Figure 5.3, shown on following page, 
we can appreciate the much smoother map of the EETL with re- 
spect to that of the ETL. This turns out to have practical advantages 
when pricing bespoke attachment and detachments or forward start- 
ing tranches. 


f(t, kj-2,kj-1) — ft, kj-1, kj) 20 => 0 
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Figure 5.3 EETL calibrated to the DJ-iTraxx 5- and 10-year tranches plus index. 


5.2 NUMERICAL RESULTS 


Our sample pre-crisis goes from 13 November 2003 to 14 June 2006 
for the CDX and from 21 June 2004 to 23 May 2006 for DJ-iTraxx. 
From Table 5.1, on next page, we note that, except for the DJ-iTraxx 
pool with a linear interpolation, in all other cases we find a solution 
where the theoretical spread exceeds the bid-ask spread by less than 
one-fifth (0.4/2) of the bid-ask range. In the case of the DJ-iTraxx 
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Table 5.1 Percentage of sample repriced outside the bid-ask range. Sample data 
for CDX range from 13 November 2003 to 14 June 2006, while for DJ-iTraxx 
from 21 June 2004 to 23 May 2006. 


CDX DJ-iTraxx 
Interpolation: Linear Spline Linear Spline 
Number of dates: 616 616 473 473 
% Misprpid-Ask >1 1.0% 2.6% 0.2% 1.3% 
% Misprpid-Ask > 1.2 0.8% 0.2% 0.2% 0.2% 
% Misprpid-Ask > 1.4 0.0% 0.0% 0.2% 0.0% 
% Misprpid-Ask > 1.6 0.0% 0.0% 0.0% 0.0% 


pool with a linear interpolation, we find only one date where all 
instruments cannot be priced within the bid-ask range: in this case 
the theoretical spread is outside the bid-ask spread by less than one- 
third (0.6/2) of the bid-ask range. 

We did not consider super-senior tranches in the pre-crisis calibra- 
tion as these tranches only became liquid enough at the end of the 
sample for our pre-crisis analysis, as in Torresetti et al. (2006a). 

In all cases we find only few dates where we could not calibrate 
an EETL surface within the bid-ask spread of each tranche plus 
the index. On the few dates where market instruments could not 
be priced exactly, the standardized mispricing is extremely small: 
the highest mispricing was obtained for the DJ-iTraxx using a linear 
interpolation: the mispricing being only 0.3 times the bid-ask spread: 
0.3 = (1.6 — 1)/2. 


5.3 SUMMARY ON EXPECTED (EQUITY) 
TRANCHE LOSS 


This methodology, first illustrated in Torresetti et al. (2006a), is 
related to Walker’s (2006) earlier work and to the formal analysis of 
the no-arbitrage properties of expected tranche loss in Livesey and 
Schlögl (2006). 
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The notion of ETL surface across maturity and tranche attachments 
is as close to a model independent notion of implied dependence as 
possible, since it focuses on one of the most direct market objects 
embedded in market quotes. 

Rather than going through implied correlation based on the arbi- 
trary assumption of a Gaussian Copula connecting defaults across 
names and leading to inconsistencies in the temporal axis, one di- 
rectly considers quantities entering the valuation formula and implies 
them from market quotes, given minimal interpolation assumptions. 
To make sure that interpolation does not interfere excessively, we 
carry out the calibration through two different interpolation tech- 
niques (linear and splines). 

While in our framework the bid-ask spreads of the instruments 
enter the target function, we aim at minimizing in order to imply the 
surface, in Walker’s (2006) framework the instruments NPVs must 
be exactly zero. By including bid-ask spreads, the no-arbitrage con- 
straints are satisfied across the vast majority of dates — in particular 
we find fewer violations of the no-arbitrage condition than in Walker 
(2006). 

The method appears to be helpful as a first model-independent 
procedure to deduce implied expected loss surfaces from market 
data, allowing one to check basic no-arbitrage constraints in the 
market quotes. 

It is also of immediate use to value tranches with non-standard 
attachments and maturities and forward starting tranches, although 
excessive extrapolation is to be avoided. 

We finally mention some further references that appeared later, 
dealing with evolutions of this technique: Parcell and Wood (2007) 
consider carefully the impact of different kinds of interpolation, 
whereas Garcia and Goossens (2007) compare the EETL between 
the Gaussian Copula and Lévy models. 


6 


A Fully Consistent Dynamical 
Model: Generalized-Poisson 
Loss Model 


The Generalized-Poisson Loss (GPL) model can be formulated as 
follows. Consider a probability space supporting a number n of inde- 
pendent Poisson processes Nj, ..., N, with time-varying, and possi- 
bly stochastic, intensities A, ..., àn under the risk-neutral measure 
Q. The risk-neutral expectation conditional on the market informa- 
tion up to time f, including the pool loss evolution up to f, is denoted 


by E,. Intensities, if stochastic, are assumed to be adapted to such 
information. 
Define the stochastic process 


Ze: = > oNj(), (6.1) 
j=l 
for positive integers &1,..., œn. In the following we refer to the 


Z, process simply as the GPL process. We will use this process as 
a driving process for the cumulated portfolio loss L,, the relevant 
quantity for our payoffs. In Brigo et al. (2006a, 2006b) the possible 
use of the GPL process as a driving tool for the default counting 
process C; is illustrated instead. 

The characteristic function of the Z, process is 


pz (u) = Bofe” ] = Bof e" | A(t)... A0) | 
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where A ;(t) := h à;(s)ds, with i = 1...n, are the cumulated in- 
tensities of each Poisson process. Now, we substitute Z, obtaining 


pz, (u) = Eo I] Bof eM | Ait)... An(t) | 
j=1 


n 


= Eo I] PN, (tA (1) UO; ) 
j=l 


which can be directly calculated since the characteristic function 
Pn; (1)|A;(t) Of each Poisson process, given its intensity, is known in 
closed form, leading to 


pz, (u) = Eo} exp | >> Aae" — 1) (6.2) 


i=l 


The marginal distribution pz, of the process Z, can be directly 
computed at any time via inverse Fourier transformation of the char- 
acteristic function of the process. The characteristic function @z,(u) 
can be explicitly calculated for some relevant choices of Poisson cu- 
mulated intensity distributions — see, for example, Brigo et al. (2006a, 
2006b). 


6.1 LOSS DYNAMICS 


The underlying GPL process Z, is non-decreasing and takes ar- 
bitrarily large values, given large enough times. The portfolio cu- 
mulated loss and the rescaled number of defaults processes are 
non-decreasing, but limited to the interval [0, 1]. Thus, we consider 
the deterministic non-decreasing function Y : N U {0} > [0, 1] and 
we define the cumulated portfolio loss process L, as 


= L 
L, := V,(Z,) := min(Z,, M^), and L, := ai (6.3) 


where 1/M’, with M’ > M > 0, is the minimum jump for the loss 
process. 
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Remark 6.1.1 Notice that the loss is bounded within the interval 
[0,1] by construction, but there is still the possibility that the loss 
jumps more than M times, where M is the number of names in 
the portfolio. If this is the case, we may check a posteriori that 
the probability of such events is negligible. This happens in all our 
examples. 


The marginal distribution of the cumulated portfolio loss process 
L, can be easily calculated. We obtain 


L, = min(Z,, M’) = Zi liz,<mM'} F M'1liz,>m') 


Since Z, has a known distribution, the distribution of L, can be easily 
derived as a by-product. The related density (defined on integer values 
since the law is discrete) is 


PLAX) = pz) lpm} + QZ: > M1 peu 


The density of L, follows directly. Also, the intensity of L;, i.e. 
the density of the absolutely continuous compensator of L, (see, for 
example, Giesecke and Goldberg, 2005), can be computed directly 
and is given by 


hy(t) = $ min(a;, (M' — Z,-))aj(0), (6.4) 
j=l 

see Brigo etal. (2006a, 2006b) for details. The intensity h goes to zero 
when Z exceeds M’, which corresponds to total loss, as expected. 
Further, if all the possible integer jump sizes between 1 and M’ are 
allowed, i.e. if; = j andn = M’, the intensity hz jumps whenever 
the cumulated portfolio loss process L jumps. The intensity jumps 
downwards, and this would seem to go in the opposite direction with 
respect to self-excitedness, which is considered a desirable feature of 
loss models in general. However, self-exciting features are embedded 
in our model, and they are embedded in the possibility of having 
several defaults at small intervals, contrary to most approaches to loss 
modelling. Consider, for example, just two names: instead of having 
the loss of one name increase the likelihood of default (intensity) of a 
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second name, we have both names defaulting together immediately. 
This embeds self-excitedness, although in an extreme way. Finally, 
to view the effect of single names on each other, one needs a more 
sophisticated formulation, based on the common Poisson shocks 
framework, leading to the GPCL model, which is analysed in Brigo 
et al. (2007). 


6.2 MODEL LIMITS 


The GPL model we introduced so far can be viewed as a particularly 
simple parameterization of the market implied loss distribution dy- 
namics. A positive feature is that the loss changes only by positive 
jumps, as should happen in any sensible loss model. Furthermore, this 
choice allows us to achieve a good calibration to market data, as we 
can see in the following section. However, we are not making explicit 
assumptions on two important issues. First, we have not addressed 
possible ways to make our model consistent with single-name dy- 
namics. Second, we have not explained how to choose a full-featured 
pool spread and recovery dynamics: a possibility would be to make 
the intensities in the Poisson processes driving Z stochastic and to 
consider more general transformations of Z to obtain the loss process. 

Since, in this book, we focus on the calibration of CDO tranches, 
which only depend on the loss marginal distribution, we avoid a 
discussion of such problems and direct the readers to Brigo et al. 
(2006a, 2006b) for an extensive analysis of candidate spread and 
recovery dynamics, and to Brigo et al. (2007) for further discussion, 
including consistency with single-names data. We point out that 
significant progress and testing in loss modelling will be possible 
only when more liquid market quotes for tranche options and forward 
start tranches are available. 


6.3 MODEL CALIBRATION 


We work with the basic GPL model specification given by the driving 
GPL process Z in (6.1), which we use in modelling the pool loss 
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through (6.3). In this basic formulation each Poisson mode N; has a 
deterministic piecewise-constant intensity A ;(t). 

Given that we have modelled the pool loss L; directly, we do not 
characterize completely the rescaled default counting process C,, but 
we give only its expectations. Indeed, such expectations are the only 
information on default counting that are implicit in index market 
quotes (2.1), whereas tranche quotes (2.4) depend only on the loss 
and not explicitly on default counting. We thus assume 


Eo[ C, | := ——E,[L,] with 0<R<1-E)[Lz,] 
(6.5) 


where the range of definition of the constant R is taken in order to 
ensure that at each time t the expected value of the rescaled number of 
defaults is greater, or equal to, the cumulated portfolio loss, and that 
both are smaller or equal to 1. Notice that we do not wish to introduce 
an explicit dynamics for the recovery rate—see Brigo et al. (2006a, 
2006b and 2007) for an initial discussion on recovery dynamics. Our 
R here can be interpreted as a sort of average recovery rate. 


6.4 DETAILED CALIBRATION PROCEDURE 


The model parameters found by the calibration procedure are the 
amplitudes a; € {m e N:m < M’} with j = 1...n, and the cu- 
mulated intensities A ;(T), which are real non-decreasing piecewise 
linear functions in the tranche maturity. 

The optimal values for the amplitudes a are selected in the follow- 
ing way: 


1. Fix the minimum jump size to 1/M’ by choosing the integer 
M'>M>O0. 

2. Find the best integer value for a, by calibrating the cumulated 
intensity A, for each value of a in the range [1, M’], all other 
modes being set to zero. 
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3. Add the amplitude œ and find its best integer value by calibrating 
the cumulated intensities A; and Ag, starting from the previous 
value for A; as an estimate, for each value of a in the range 
[1, M°], 

4. Repeat the previous step for œ; with i = 3, and so on, by cali- 
brating the cumulated intensities A,,..., A;, starting from the 
previously found A;,..., A;—; as the initial estimate, until the 
calibration error is under a given threshold, or until the intensity 
A; can be considered negligible. 

5. Check a posteriori that the probability of having more than M 
jumps is negligible and that the value of is within the arbitrage- 
free range given in (6.5). 


The objective function f to be minimized in the calibration is the 
squared sum of the errors shown by the model to recover the tranche 
and index market quotes weighted by market bid-ask ai 


A 
f(a, A) = Lei < c= -10a = (6.6) 


where the x;, with i running over the market quote set, are the index 
values So for DJ-iTraxx index quotes, and either the index periodic 
premiums S: 2 or the upfront premium rates U ^? for the DJ-iTraxx 
tranche quotes. 


6.5 CALIBRATION RESULTS 


The GPL model is calibrated to the market quotes observed weekly 
from 6 May 2005 to 18 October 2005. We take R = 30%, following 
Albanese et al. (2006), as the reference value for the recovery rate 
in the DJ-iTraxx Europe market for spot and forward contracts. The 
quality of our calibration below is not altered if we select a value R = 
40% resembling the recovery typically used in simplified quoting 
mechanisms in the market- see Brigo et al. (2006a, 2006b, 2007) 
for some examples. We start with M’ = 200, corresponding to a 
minimum loss jump size of 50 bp. 
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Table 6.1 DJ-iTraxx index and tranche quotes in basis points on 13 May 2005, 
along with the bid-ask spreads. Index and tranches are quoted through the periodic 
premium, whereas the equity tranche is quoted as an upfront premium (See 
Chapter 2 for the market quoting conventions.) 


Maturities 
Att—Det 3 yr 5 yr Tyr 10 yr 
Index 38(4) 54(1) 65(3) 77(2) 
Tranche 0-3 2060(100) 4262(118) 5421(384) 6489(124) 
3-6 72(10) 173(68) 398(40) 590(20) 
6-9 28(6) 57(6) 141(17) 188(15) 
9-12 13(2) 31(5) 72(20) 87(15) 
12-22 3(1) 21(3) 42(13) 60(10) 


Consider, as a first example, the calibration date 13 May 2005 (see 
Table 6.1). In Table 6.2 we list the calibration result and the values of 
the calibrated parameters. The calibration errors € are very low for 
all maturities. Notice that a calibration error smaller than 1 means 
that the difference between the market quote and the model price is 
smaller than the bid-ask spread. 


Table 6.2 Left-hand side: calibration error €; in (6.6) (calculated with respect to 
the bid-ask spread) for tranches quoted on 13 May 2005. Right-hand side: 
cumulated intensities, integrated up to tranche maturities, of the GPL model with 
M' = 200. Each row corresponds to a different Poisson component with jump 
amplitude a. Recovery: 30% 


Maturities A(T) 


Att-Det 3yr Syr 7yr 10yr a 3yr Syr 7yr 10yr 


1 1.955 3.726 4.464 7.694 
3 0.000 0.062 0.305 0.305 
Tranche 0-3 0.0 O1 0.2 —0.2 8 0.016 0.033 0.109 0.109 
3-6 0.0 0.0 —0.2 0.0 12 0.004 0.013 0.026 0.026 
6-9 0.0 00 —0.3 0.1 19 0.006 0.006 0.017 0.017 
9-12 —0.1 0.1 —0.1 04 72 0.000 0.009 0.026 0.049 
12-22 0.0 0.0 —0.2 —0.3 185 0.000 0.002 0.002 0.008 


Index 0.0 —0.1 03 0.0 
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Table 6.3 Same as Table 6.1 on 11 October 2005 


Maturities 

Att—Det 3 yr 5 yr Tyr 10 yr 

Index 23(2) 38(1) 47(1) 58(1) 
Tranche 0-3 762(26) 3137(26) 4862(76) 5862(74) 
3-6 20(10) 95(1) 200(3) 515(10) 

6-9 7(6) 28(1) 43(2) 100(4) 

9-12 12(2) 27(4) 54(5) 

12-22 71) 13(2) 23(3) 


Consider, as a second example, the calibration date October 11, 
2005 (see Table 6.3). We list in Table 6.4 the calibration results and 
the values of the calibrated parameters. The calibration errors show 
that the 10-year equity tranche is not correctly priced. We find such 
mispricing in many calibration examples, in particular after October 
2005. 

We then try as minimum loss jump size 1/M’, besides 50 bp, the 
values 2 bp and 10 bp, corresponding, respectively, to M’ equal to 
5000 and 1000. Table 6.5, on next page, shows that the 10-year ma- 
turity tranches (the most difficult to calibrate in our experience) are 
stable through the three different loss sizes, suggesting that going 


Table 6.4 Same as Table 6.2 on 11 October 2005 (the 3-year maturity quotes 
lack two tranches) 


Maturities A(T) 


Att-Det 3yr Syr 7yr 10yr a 3yr 5yr 7yr 10yr 


Index 0.0 00 01l 01 1 0.441 2.498 4.466 7.555 
2 0.435 0.435 0.435 0.671 
Tranche 0-3 —0.1 0.1 —1.2 21 11 0.004 0.023 0.023 0.023 
3-6 —0.1 —0.1 0.3 —1.0 22 0.000 0.001 0.006 0.030 
6-9 0.0 —0.1 03 09 29 0.000 0.000 0.001 0.001 
9-12 0.4 —0.8 —0.8 32 0.001 0.004 0.004 0.004 
19299 0.0 0.0 0.0 192 0.000 0.001 0.005 0.011 
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Table 6.5 GPL calibration error for different minimum loss sizes 1/ M’ 
with respect to the bid-ask spread for 10-year tranches on 11 October 
2005. Recovery: 30% 


Att—Det 50 bp 10 bp 2 bp 

Tranche 0-3 2.1 1.8 1.8 
3-6 —1.0 —1.0 —1.0 

10-yr maturity 6-9 0.9 0.9 0.9 
9-12 —0.8 —0.9 —0.8 

12-22 0.0 0.2 0.0 


below 50 bps does not add much flexibility to the model. This is 
confirmed by further tests, and in particular the difference between 
the M’ = 1000 and the M’ = 5000 calibrations is always small. Fur- 
thermore, the behaviour of the mean calibration error, i.e. of the mean 
of the absolute values of the €; across time and quotes, for the three 
different choices of M’, is quite similar and within 1. 

We noticed also that, as the minimum jump size decreases 
(“granularity” increases), the loss distribution becomes noisier, due 
to the presence of small amplitudes. Furthermore, very small modes, 
appearing when the minimum jump size is as small as a few basis 
points, may violate the requirement that the loss process jumps less 
than M times (see Remark 6.1.1). We also tried calibrations with M’ 
less than 200, i.e. with minimum loss jump greater than 50 bp. In 
this case the calibration error grows quickly. Indeed, the minimum 
jump size, in this case, becomes greater than the typical portfolio loss 
given when one name defaults. 50 bp seems then to be a reasonable 
reference value. 

Also, the values of the Poisson amplitudes are quite stable across 
the calibration dates. Indeed, in six months we observe at most four 
changes in their values, as shown in Table 6.6, on following page. 

The loss distribution implied by the GPL model is multimodal and 
the probability mass moves towards larger loss values as the maturity 
increases. These features are shared by different approaches. For 
instance, static models, such as the implied default-rate distribution 
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Table 6.6 Values of the Poisson amplitudes œ/ M’, for different values of the 
minimum loss jump 1/M’. Only the calibration dates between 6 May 2005 and 18 
October 2005 where the a/M’ values change are listed 


50 bp Poisson’s amplitudes 


Date 1 2 3 4 5 6 7 


06-May-05 0.50% 1.50% 4.00% 6.00% 9.50% 39.50% 92.50% 
02-Sep-05 0.50% 1.00% 4.00% 5.50% 12.50% 39.00% 100.00% 
11-Oct-05 0.50% 1.00% 5.50% 11.00% 14.50% 16.00% 96.00% 


10 bp Poisson’s amplitudes 


Date 1 2 3 4 5 6 7 


06-May-05 0.10% 1.50% 4.60% 5.90% 9.60% 39.60% 53.00% 
05-Aug-05 0.20% 1.10% 1.40% 8.10% 11.30% 49.00% 62.40% 
11-Oct-05 0.10% 0.70% 1.00% 6.30% 11.50% 14.50% 93.70% 


2 bp Poisson’s amplitudes 


Date 1 2 3 4 5 6 7 


06-May-05 0.02% 1.50% 5.26% 9.64% 17.58% 39.64% 99.78% 
12-Aug-05 0.38% 1.06% 1.14% 7.38% 12.24% 41.34% 99.80% 
03-Oct-05 0.02% 0.98% 1.16% 7.52% 9.74% 43.34% 65.16% 
11-Oct-05 0.16% 0.68% 1.00% 6.30% 10.98% 14.46% 94.90% 


in Torresetti et al. (2006c), suggest multimodal loss distributions, as 
we have seen also here in our earlier chapters on the Implied Copula. 
The evolution of the implied loss distribution is shown in Figure 6.1 
on next page. 

The dynamic credit correlation model by Albanese et al. (2006) 
shows implied loss distributions whose modes tend to group as the 
maturity increases, leading to a distribution approaching normality. 
The GPL model reproduces this behaviour (see, for example, Brigo 
et al., 2006a, 2006b). 
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Figure 6.1 Loss distribution evolution of the GPL model with minimum jump 
size of 50 bp at all the quoted maturities up to 10 years, drawn as a continuous line. 


7 
Application to More Recent Data 


and the Crisis 


In this chapter we check whether the critical features we discussed 
about implied correlation and the subsequent elements coming from 
more advanced models are still present in-crisis, after mid-2007. We 
will observe that the features are still present and often amplified in 
the market after the beginning of the crisis. 


7.1 COMPOUND CORRELATION IN-CRISIS 


We begin by extending the sample of the compound correlation 
invertibility analysis to see if the non-invertibility arose again during 
the market turmoil in 2007-2009 (see Figures 7.1 to 7.4 on pages 88 
to 91). As before, we resort to the homogeneous pool model. 

In Figures 7.2 and 7.4 we notice how, since the credit crunch 
began in Summer 2007, the problematic tranches changed. In fact, 
the non-invertibility affected 


(1) For the DJ-iTraxx, 10-year 12—22% tranche (see Figure 7.2); 
(2) For the CDX pool the 10-year 10-15% and more marginally the 
10-year 15-30% tranche (see Figure 7.4). 


Another relevant issue, besides the non-invertibility of compound 
correlation, is its non-uniqueness when invertible. We have seen 
in Figure 3.1 that on 3 August 2005 for some tranches the map 
between correlation and theoretical tranche spread is not mono- 
tonic. In particular, we have seen how this resulted for the DJ-iTraxx 
10-year 6-9% and the DJ-iTraxx 10-year 9-12% tranches. 
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Figure 7.1 DJ-iTraxx 5-year compound correlation invertibility. Market spread 
(black line) versus minimum and maximum tranche spread obtained varying the 
compound correlation between 0 and 1 (grey lines). Dots are highlighting the dates 
and market spread that were not invertible (black line lies outside the area within 
the two grey lines). 
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Figure 7.2 DJ-iTraxx 10-year compound correlation invertibility. Same as Fig- 


ure 7.1. 
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Figure 7.3 CDX 5-year compound correlation invertibility. Same as Figure 7.1. 
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Figure 7.4 CDX 10-year compound correlation invertibility. Same as Figure 7.1. 
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Despite the map for these two tranches not being monotonic, on 
this particular date we did not encounter the problem of having more 
than one compound correlation that implied the same theoretical 
spread. 

In Figures 7.5 to 7.8, on pages 93 to 96, we check whether 
there have been dates where a given tranche market spread could 
be mapped into more than one compound correlation. We plot the 
market spread and highlight with dots the dates when the market 
spread was invertible into two different compound correlations. 

We note that the problem of the non-uniqueness of the compound 
correlation is particularly acute for several mezzanine tranches for 
both maturities (5 and 10 years), and for both regions (DJ-iTraxx 
and CDX). Also, for a given tranche the problem of non-uniqueness 
of the compound correlation can span a relevant time window: for 
example, the 5-year DJ-iTraxx 3—6% compound correlation has not 
been uniquely invertible from March 2005 to October 2007. 

At this point we aim at addressing the question of the economic 
impact of the non-uniqueness of the compound correlation. 

Let us take, for example, the DJ-iTraxx 5-year 6-9% tranche on 
2 July 2008. In Figure 7.9 on page 97 we see that this is a date where 
the tranche market spread of 261 bp was not uniquely invertible: the 
same market spread could be obtained with more than one compound 
correlation. 

In Figure 7.9 we plot the market spread (flat grey line) and the 
theoretical spread (dotted black line) as a function of the compound 
correlation. The same tranche spread could be obtained with a com- 
pound correlation of either 2.41% or 95.60%. 

In Figure 7.10 on page 97 we plot the ETL as a function of time 
for the two correlations, 2.41% or 95.60%. We note the completely 
different profile of the ETL, which in turn translates into sensibly 
different forward start spreads and zero coupon prices for the same 
DJ-iTraxx 5-year 6-9% tranche. 

The payoff of the zero coupon tranche is simply the outstanding 
notional left of the tranche at maturity. Thus, by assuming deter- 
ministic default-free interest rates, the price at time 0 of the zero 
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Figure 7.5 DJ-iTraxx 5-year compound correlation uniqueness. Black dots are 
highlighting the dates where more than one compound correlation could reprice the 
tranche market spread. 
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Figure 7.6 DJ-iTraxx 10-year compound correlation uniqueness. Same as Fig- 


ure 7.5. 
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Figure 7.7 CDX 5-year compound correlation uniqueness. Same as Figure 7.5. 
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Figure 7.8 CDX 10-year compound correlation uniqueness. Same as Figure 7.5. 
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Figure 7.9 Market spread (flat line) and theoretical spread as a function of the 
compound correlation (dotted line) of the DJ-iTraxx 5-year 6-9% tranche on 2 July 
2008. 
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Figure 7.10 Expected loss of the DJ-iTraxx 5-year 6-9% tranche for a compound 
correlation of 95.60% (continuous line) and for a compound correlation of 2.41% 
(dotted line), both consistent with the market spread of the DJ-iTraxx 5-year 6-9% 
tranche on 2 July 2008. 
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As for the forward start tranche, the forward tranche protection 
buyer will begin paying a premium to the protection seller from the 
forward start date T, to the tranche maturity T = Tp. The outstanding 
notional on which these payments are calculated will be the average 
outstanding notional of the tranche on each interest rate period. The 
protection seller will, instead, pay the protection buyer all default 
payment accrued by the tranche from T, to Tp. 

This means that if all the tranche notional is wiped out by the 
forward start date — i.e. the portfolio loss by that date has reached the 
detachment — no flows are exchanged between the protection buyer 
and seller. Thus, the fair spread of a forward start tranche will be 
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In Table 7.1 we compare the default leg price, the zero coupon tranche 
price and the forward start (20 June 2009) tranche spread associated 
with the two compound correlations. While both correlations set the 
tranche premium leg value equal to the tranche default leg value, the 


Table 7.1 Default leg price, zero coupon price and forward start (20 June 2009) 
spread of the 5-year DJ-iTraxx 6-9% tranche associated to the two compound 
correlations 95.60% and 2.41% consistent with the market tranche spread on 2 
July 2008 


Compound Default Zero coupon Forward start 
correlation leg price tranche price tranche spread 
2.41% 11.23% 10.71% 314 bp 
95.60% 10.80% 9.43% 258 bp 
Absolute difference 0.43% 1.28% 56 bp 


Relative difference 3.98% 13.57% 21.71% 
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tranche riskiness, as measured by the value of the default leg, need 
not be equal for both correlations. 

From Table 7.1 we also note that while the default leg price dif- 
ference is only 43 bp, or 3.95% in relative terms (11.23%/10.80% 
—1), the differences in relative terms for the zero coupon tranche 
price and the forward start tranche spread are more significant: re- 
spectively 13.53% and 21.7%. 

Thus besides being difficult to use for bespoke tranches, compound 
correlation is often non-invertible or admits multiple solutions. When 
multiple solutions are present, they may imply different values for 
zero coupon tranche prices and forward start tranche spreads (and 
these are only two examples). 


7.2 BASE CORRELATION IN-CRISIS 


In section 3.7 we pointed to the solutions that the base correlation 
provided for some of the limitations of compound correlation. We 
also pointed out the particular care one needs to take when dealing 
with base correlation given its intrinsic inconsistency, using two dif- 
ferent distributions to calculate the same tranche value. In particular, 
as a result of this inconsistency, we have observed negative expected 
tranche losses coming from base correlation. As before, we resort to 
the homogeneous pool model. This renders the analysis a little less 
realistic, since the heterogeneous model is used in practice. However, 
difficulties in calibration of base correlation with the homogeneous 
model give a very strong alarming signal per se, and the heteroge- 
neous version in such a case is not expected to perform any better. 

Overall, the base correlation surpassed the compound correlation 
and prevailed over more sophisticated dynamic loss models pre-crisis 
because: 


e it could be calibrated to almost all market conditions experienced 
pre-crisis; 

e it resulted in a smooth map, thus allowing for pricing of bespoke 
tranches via interpolation; 
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it could be used to define sketchy correlation mapping methods 
used for pricing tranches on bespoke CDO portfolios, where no 
other correlations were available — see, for example, Bear Stearns 
and Lehman Brothers’ reports by Reyfman et al. (2004) and Baheti 
and Morgan (2007); 

e in its heterogeneous version it provided traders with easy to calcu- 
late hedge ratios. 


In-crisis this perception was quite shaken. For example, the fi- 
nite pool heterogeneous one-factor Gaussian Copula base correlation 
could often not be calibrated to the 5-year CDX Series 9 tranches 
during 2008.! 

Let us keep in mind that the base correlation associated with the de- 
tachment B is the correlation for an equity tranche with a detachment 
equal to B. From equation (2.3) we note that the expected loss of the 
tranche 0, B, in turn, enters with a positive coefficient the expected 
loss of the tranche A, B. From our example in Figure 3.7, we see that 
an increasing 0, B detachment base correlation, when keeping fixed 
the attachment 0, A base correlation, results in a decreasing tranche 
spread. 

When the 15-30% tranche spread cannot be calibrated, this hap- 
pens because the 30% detachment base correlation cannot be in- 
creased enough: the minimum attainable theoretical 15-30% tranche 
spread obtained for a 30% detachment base correlation of 100% is 
above the 15—30% market tranche spread. This, in turn, occurs be- 
cause the 15% detachment base correlation is already so high that 
the 100% upper bound for correlation values becomes a binding 
constraint. 

In this respect we plot, in Figure 7.11 on following page, the CDX 
5-year Series 9 15-30% tranche market spread and the minimum 
tranche spread that could be calibrated using the finite pool heteroge- 
neous one-factor Gaussian Copula base correlation from September 
2007 to May 2009. 


' We have limited our attention to the Series 9 as the liquidity on the subsequent series has been 
somewhat reduced. 
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Figure 7.11 Market spread and minimum attainable spread via heterogeneous 
base correlation calibration for the 5-year CDX 15-30% tranche. 


We have thus seen that the heterogeneous base correlation had 
calibration problems on those dates where the 100% upper bound 
for the correlation became a binding constraint. If we move to a 
homogeneous assumption, the homogeneous copula could calibrate 
all the capital structure where the heterogeneous copula could not. 

In fact, allowing heterogeneity may change considerably the shape 
of the loss distributions, as can be seen in Figure 7.12 on next page 
where we compare the 5-year CDX pool loss distribution implied 
by the same correlation parameter for the homogeneous and het- 
erogeneous copulas. The correlation parameter used is the one that 
calibrates the 0-3% equity tranche spread on 7 December 2007 under 
a heterogeneous model that is consistent with the constituents CDS 
quotes. 

In Figure 7.13 on next page we plot the 5-year CDX homogeneous 
and heterogeneous base correlations calibrated to the market spread 
on 7 December 2007. We note that the heterogeneous base correlation 
could not calibrate the 15-30% tranche and we also note that the 
homogeneous base correlation is smaller than the heterogeneous 
base, so that the 100% upper bound for correlation was not reached. 
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Figure 7.12 Pool loss distribution computed via a finite pool heterogeneous 
(black line) and homogeneous pool (grey line) one-factor Gaussian Copula. 
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Figure 7.13 Heterogeneous and homogeneous base correlation calibrated to the 
5-year CDX tranche spreads on 7 December 2007. 
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Summarizing, the heterogeneous version of the base correlation 
failed during the recent credit crisis as market participants have not 
been able to calibrate it for a considerable amount of time. The 
homogeneous version of the base correlation instead has been able 
to calibrate both 5- and 10-year maturities and both DJ-iTraxx and 
CDX indices throughout the crisis. 

This does not mean that the latter has been immune from problems. 
We are going to highlight those problems in order to make a fair 
assessment of the homogeneous pool base correlation versus the 
Implied Copula and the GPL in-crisis. 

Most problems in base correlation calibration concerned senior or 
super-senior tranches. Within the base correlation framework, the 
super-senior tranche pricing is a by-product of the calibration to 
the rest of the capital structure and the index. For example when 
calibrating the 22-100% tranche value, it will be necessary, at the 
end of the iterations through the capital structure dictated by the 
base correlation calibration procedure, to compute the following 
expectation: 


ni 7 1 : . 
[Jona -rapele 


-0.22 E | 27°? | (6022)] 


where the parameter to move in order to calibrate the tranche value 
is p1. 

The problem is that the expectation of the entire pool loss is obvi- 
ously not sensitive to correlation, while being related to index quotes. 
It is in this sense that we would obtain the super-senior tranche spread 
as a by-product of the base correlation calibration procedure and the 
index. 

We define the relative mispricing as the amount by which the 
theoretical spread falls outside the bid-ask spread divided by the 
mid-market spread. In Figure 7.14, on following page, we plot the 
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Figure 7.14 Relative mispricing of the CDX 5-year 30-100% when the rest of 
the capital structure is calibrated with a finite homogeneous pool base correlation 
(when the 30% base could not be calibrated the mispricing was set to zero). 


relative mispricing, the same metric that will be used to judge the 
goodness of the GPL calibrations in-crisis, for the CDX 5-year 30- 
100% super-senior tranche. We notice that the mispricing has been 
quite relevant: hovering in the 10-20% range most of the time with 
spikes in the 50-60% range. 

Finally, we mention that the deterministic recovery assumption, 
while being computationally very convenient, does not help in the 
above setup. This has been addressed in the base correlation frame- 
work by Amraoui and Hitier (2008) and Krekel (2008), among oth- 
ers, who introduce random recovery. However, even with this update, 
base correlation remains flawed and may still lead to negative loss 
distributions, and anecdotally there are dates in-crisis where the het- 
erogeneous pool base correlation with random recovery cannot fit the 
market if the said random recovery is imposed to be consistent with 
a single-name CDS valuation. For a formulation of random recovery 
as a function of the systemic factor in the Impiled Copula framework, 
see Remark 4.1.1. 
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7.3 IMPLIED COPULA IN-CRISIS 


Here we analyse the Implied Copula framework described in Chap- 
ter 4, showing the implied loss distribution across time, in order to 
include the earlier examples as special cases of our analysis of the 
implied loss through history, both pre- and in-crisis. 

In Figure 7.15 we show the implied distribution for the default 
counting process, calibrated to the CDX and DJ-iTraxx tranches (the 
entire capital structure except the super-senior tranche) and the index 
for both the 5- and 10-year maturities from March 2005 to January 
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Figure 7.15 Implied distribution calibrated with the Implied Copula, via a two- 
stage optimization as in Torresetti et al. (2006c), to the CDX and DJ-iTraxx tranches 
for both the 5- and 10-year maturities calibrated from March 2005 to January 2009. 
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2009. We note how, since the beginning of the crisis, the probability 
mass underlying the implied loss distribution shifted towards a higher 
number of defaults for both the DJ-iTraxx and the CDX and for the 
5- and 10-year maturities. This can be interpreted as an increased 
perceived default riskiness in general. 

Also, we note how, since the start of the crisis, the probability mass 
associated with a catastrophic or Armageddon event — i.e. the default 
of the entire pool of credit references — has increased dramatically. 
This is particularly visible for the 10-year default rate distributions for 
both DJ-iTraxx and CDX. The increased probability of Armageddon 
has also been pointed out in Morini and Brigo (2007, 2009) in the 
context of credit index options. 

In Figures 7.16 and 7.17 on pages 107 and 108 we show a feature of 
the implied distribution that cannot be immediately appreciated when 
inspecting Figure 7.15: the stability of the location of the modes in 
the implied loss distribution through time. Taking a snapshot of the 
implied distribution calibrated to the CDX 5-year tranches approx- 
imately every five months, we note how the mode associated with 
about 20 names is quite constant in its location across all calibra- 
tion dates. In particular, on recent data, we note how even although 
the distribution shifts to the right, thus pricing an increased risk of 
default of the pool, the location of the mode does not seem to be 
affected. In particular, we examine data on March 2008 and Jan- 
uary 2009: despite the massive shift to the right of the distribution, 
the non-parametric optimization we implemented still shows a clear 
mode that would be associated with a cluster of defaults amounting 
to about 20 entities. This observation leads almost naturally to the 
GPL dynamic framework outlined in section 7.5. 


74 EXPECTED TRANCHE LOSS 
SURFACE IN-CRISIS 


The performance of the Expected Tranche Loss (ETL) surface model 
in-crisis was not as egregious as pre-crisis. This was in part due to the 
fact that with the super-senior tranche available, the system became 
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Figure 7.16 Implied distribution calibrated with the Implied Copula, via a two- 
stage optimization as in Torresetti et al. (2006c), to the CDX 5-year on-the-run 
tranches from June 2005 to June 2007. See Figure 7.17 for data after June 2007. 


overdetermined. In fact the capital structure was completed? and the 
super-senior expected loss did not only have to fit the index spread, 


2 Since the end of 2005 the liquidity of the super-senior tranche allowed us to include it in the 


set of calibration instruments. 
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Figure 7.17 Implied distribution calibrated with the Implied Copula, via a two- 
stage optimization as in Torresetti et al. (2006), to the CDX 5-year on-the-run 
tranches from October 2007 to January 2009. See Figure 7.16 for data before 
October 2007. 


via the deterministic recovery, but had to retrieve both the index and 
the super-senior tranche. 

Adding the super-senior tranche to the set of market instruments 
in the optimization of (5.4) does not add any degrees of freedom in 
the variables whose function is being optimized, while adding two 
instruments: one for each maturity. 

With the addition of the super-senior tranche the unknown ETL of 
the 22-100% tranche, for example, will not be calibrated to the index 
assuming a deterministic recovery, but will be calibrated instead 
directly to the super-senior tranche spread. The index will then be a 
by-product of the calibration of all other tranches if one assumes a 
deterministic recovery. 
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This procedure could well result in a theoretical tranche-derived 
index spread that lies outside of the actually quoted index bid-ask. 
This is why, in the joint optimization of (5.4) to all tranches and the 
index, we will seek ETLs such that all tranches and the index are 
priced within the bid-ask spreads for both maturities. 

In Figure 7.18 we report the relative mispricing of the 5- and 10- 
year tranches and indices for the DJ-iTraxx Series 9 and the CDX 
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Figure 7.18 Maximum absolute standardized mispricing, i.e. absolute value 
of (7.3), across the 5- and 10-year tranches and indices resulting from the cali- 
bration on each date of the ETL surface. 
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Series 10 in-crisis when, in the optimization of (5.4), we are including 
the super-senior tranche. We note that the quality of the calibration is 
somewhat reduced: pre-crisis only 0.2% and 0.8% of all calibration 
dates for, respectively, the DJ-iTraxx and the CDX resulted in a 
maximum standardized mispricing across all instruments that was 
greater than 20% of half the bid-ask spread when interpolating the 
ETLs linearly to compute the tranche and index NPVs: see Table 5.1. 
In-crisis this figure increases to 1.4% and 1% of all calibration dates 
for, respectively, the DJ-iTraxx and the CDX: considering all that we 
went through in the last couple of years this is quite a remarkable 
achievement, even though it is a sensible deterioration with respect 
to pre-crisis calibrations. 

In all cases but one, the largest mispricing arose from the 5-year 
super-senior tranches and the index. The only exception was on 22 
January 2009, where the largest mispricing arose from the 10-year 
super-senior tranche and the index. In Figure 7.19 on following page 
we see that if the calibration results in a standardized mispricing 
greater than zero, then the mispricing is positive for the index and 
negative for the super-senior tranche.* 

We limit our attention to the DJ-iTraxx Series 9 and the CDX 
Series 10 because, during the crisis, the liquidity on the CDO market 
completely disappeared for quite an extended period on the on-the- 
run series.* Given the somewhat reduced liquidity of the CDO market 
in-crisis, we will only calibrate tranches and indices with the most 
liquid maturities: i.e. 5 and 10 years. 

A positive standardized mispricing for the index means that the 
theoretical index spread is above the offer. Conversely, for the super- 
senior tranche a negative standardized mispricing means that the 
theoretical super-senior tranche spread lies below the bid. 


3 We do not show all other tranche mispricings as they have the same sign of the super-senior 
and are usually much smaller in absolute value. 

4 The interest of market participants for the DJ-iTraxx Series 9 and the CDX Series 10 never 
disappeared because of the need to hedge legacy CDO portfolios whose rating composition was 
more in line with the off-the-run series rather than the on-the-run series. Indeed, in the on-the-run 
series, downgraded names are excluded when rebalancing the index. 
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7.4.1 Deterministic piecewise constant recovery rates 


Below we show a way to correct the overdetermined feature of our op- 
timization in (5.4), where, when adding the super-senior tranches, we 
seek 12 unknowns to calibrate 14 market instruments. This overde- 
termination resulted in the problems illustrated in Figure 7.19. 
Intuitively, we will try to increase the theoretical calibrated—-ETL 
index spread so that it falls within the quoted index bid-ask spread. 
Given the fact that the index default leg NPV is driven by the prices 
of the tranches spanning the entire capital structure so that we cannot 
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Figure 7.19 Standardized mispricing (7.3) of the 5-year super-senior tranche and 
index resulting from the calibration on each date of the ETL surface. 
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change it, we will focus instead on the index premium leg and, in 
particular, on the index Dyo. 

To this extent we note that the outstanding notional of the index, 
entering the Dyo;, is reduced by the pool default rate and not by 
the pool loss rate. In other terms, the recovery rate does not enter 
the Dvo; of the index, while entering its default leg (this is different 
from the tranches, where the recovery rate also enters the tranche 
Dyo). In order to increase the index Dyo; and thus, given the same 
default leg NPV, reduce the theoretical index spread so that it falls 
below the offer, we will introduce a deterministic piecewise constant 
recovery rate. Our modification to the ETL approach is to calibrate all 
tranches, including the super-senior, thus spanning the entire capital 
structure, plus the index, by implying a deterministic recovery that is 
piecewise constant in time over the grid {To, T5yr, Tioyr} (two steps). 
We set 


ad E R,)dC, =dL, > Ri = Rsyr Yt as T5yr , R; = Rioyr Vt > T5yr 


with Rsyr and Rioyr deterministic constants to be determined. 
Furthermore, the objective function will be changed to penalize 
deviations of the stepwise constant calibrated recovery from 40%, the 
market standard when pricing the recovery of senior unsecured CDSs. 
The trade-off factor 100 in (7.1) is intended to move the optimization 
towards small mispricings: i.e. privileging small mispricings over 
small departures of recovery from 40%.° Finally, we notice that 
the mispricing we define below in (7.3) is different from the one 
introduced pre-crisis in (5.6). This is done to give more visibility to 
the impact of bid-ask price information. While in our definition pre- 
crisis we would be within the bid-ask window when the mispricing 
measure (5.6) was smaller than 1 in absolute value, with our new 
measure (7.3) the calibration is within the bid-ask interval if the 
mispricing is zero. This way we immediately notice that when the 


5 A mispricing of 0.01 times half the bid-ask spread will increase the objective function by 
0.01 = 0.01? x 100, whereas a deviation of 0.01 of the implied recovery from 40% will increase 
the objective function by 0.0001 = 0.017. Thus mispricing any instrument by 1% the bid-ask 
spread will be 100 times more penalizing than moving the implied recovery of 1% away from 40%. 
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model calibrates outside the bid-ask spread, the mispricing becomes 
non-zero. 


argmin 100 y >, (MiserS 1024’? 


{R5yr, Rioyr, f Syr,3%) seers fAOyr, 100%)} T A,B 


+ (Rr-04" 01) 
T 


subject to: 
0< ft,k)<1l 
fti, k) = f(ti-1,k) (7.2) 
f(t, kj-1) = f(t, kj) 

where 


A, B,th A, B,bid A,B,ba : A, B,th A, B,bid 
(Sé theor Sé ji )/S6 a if Sé theor < Sé j 
A,B _ ss fn Pe 
MosprS1pz = (Sone _ Sener if Spal = Saree 
0 otherwise 


(7.3) 


Further, particularly on all dates where the optimization set out 
in (5.4) results in a standardized mispricing of any tranche or in- 
dex for any maturity larger than zero, i.e. the theoretical spread of 
any of the tranches or index for any of the two considered maturities 
is outside the bid-ask spread, we will run the optimization set out 
in (7.1) where, with respect to the optimization (5.4), we will intro- 
duce the piecewise constant recovery rate, flat between time 0 and the 
5-year maturity and then between the 5- and 10-year maturities. This 
will give us two additional degrees of freedom in the optimization. 

More specifically, for each market date we run a two-stage opti- 
mization where the second stage with more parameters (recoveries) 
is entered only when any of the instruments is not priced within the 


6 As already done for the optimization set out in (5.4), in the optimization set out in (7.1) we will 
multiply the standardized mispricing of the index by 2 in order to facilitate a much finer repricing 
of the index. 
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Figure 7.20 Stepwise constant recovery rates implied from tranches and index 
market quotes via the optimization set out in equation (7.1). 


bid-ask spread in the first stage. This way we could calibrate all dates 
without any mispricing. 

In Figure 7.20 we plot, for each date, the calibrated piecewise 
constant recovery: the 5- and 10-year implied recoveries. We note 
that on those dates where a single fixed recovery was not able to 
calibrate all tranches and the index, the above calibration has been 
implying a lower recovery than 40% for the first 5 years. This is in 
line with the evidence that recovery rates depend on the business 
cycle, recovery rates being lower during periods of recession. It is 
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also consistent with the forecasts of a slower than expected economic 
recovery as we move out of the current crisis. 


7.5 GENERALIZED-POISSON LOSS MODEL 
IN-CRISIS 


We resort now to the GPL model, to see how this performs in-crisis. 
We change the model formulation slightly in order to have a model 
that is more in line with the current market, while maintaining all 
the essential features of the modelling approach. Compared with the 
model described in Chapter 6, we introduce the following modifica- 
tions. We fix the jump amplitudes « a priori. To do this, we reason 
as follows. Fix the independent Poisson jump amplitudes to the lev- 
els just above each tranche detachment, when considering a 40% 
recovery. 

For the DJ-iTraxx, for example, this would be realized through 
jump amplitudes a; = a;/125, where 


125 - 0.03 125 - 0.06 
as = roundup | ————— ]}, a = roundup | ———— J, 


(1 — R) (1 — R) 
125 - 0.09 125 -0.12 
a7 = roundup (a) , œg = roundup (F) ; 
125 -0.22 
œg = roundup (a) » Ajo = 125 


and, in order to have more granularity, we add the sizes 1, 2, 3, 4: 
a, = l, a2 = 2, a3 =3, a4 = 4. 

In total we have n = 10 jump amplitudes. We then slightly modify 

the obtained sizes in order to account also for CDX attachments that 


are slightly different. Eventually we obtain the set of amplitudes 


125 - a; =a; € {1, 2,3, 4, 7, 13, 19, 25, 46, 125} 
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Given these amplitudes a, we obtain the default counting process 
fraction as 


n—-1 
Cr = Lw,y=0e + LiNo), E := min (x a Nj(t), ) : 
i=l 
Now let the random time 7 be defined as the first time where 
yo", a: N;(t) reaches or exceeds the relative pool size of 1. 


t= otf : $ ajNi(t) = i 
i=l 
We define the loss fraction as 


Li := lesa — R)I e=) + lec) 
x [C1 — R) lime + linao 1 — R Ze)] 
= l>a — Rlww=oe + lesa — Ree) (7.4) 


Notice that in this way whenever the Armageddon component N, 
jumps the first time, the default counting process Č, jumps to the 
entire pool size and no more defaults are possible. 

Furthermore, whenever the Armageddon component N, jumps the 
first time, we will assume that the recovery rate associated with 
the remaining names defaulting in that instant will be zero. The pool 
loss, however, will not always jump to | as there is the possibility that 
one or more names have already defaulted before the Armageddon 
component N,, jumped, and they defaulted with recovery R. If at a 
given instant ¢ the whole pool defaulted, i.e. C; = 1, this may have 
happened in two ways: 


e N, jumped by t. In this case the portfolio has been wiped out with 
the contribution of an Armageddon event. Notice that in this case 
dC, = dL, if N, jumps at t. In fact the recovery associated with 
the pool fraction defaulting in that instant will be equal to zero. 

e N, has not jumped by f. In this case the portfolio has been wiped 
out without the contribution of an Armageddon event but because 
of defaults of more small or big sectors that do not comprise the 
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whole pool. Notice that, in this case, the loss is less than the whole 
notional, as all these defaults had recovery R > 0. 


In this way, whenever N, jumps at a time when the pool has not 
been wiped out yet, we can rest assured that the pool loss will be 
above | — R. We do this because the market in 2008 has been quoting 
CDOs with prices assuming that the super-senior tranche would be 
impacted to a level that is impossible to reach with fixed recoveries 
at 40%. Anecdotically, there was a market for the DJ-iTraxx 5-year 
60-100% tranche on 25 March 2008, quoting a running spread of 
24 bp bid. 

We know how to calculate the distribution of both C, and L,, given 
that: 


e the distribution of ¢, = min paa ai N;(t), 1) is obtained running 
a “reduced" GPL, i.e. a GPL where the jump N, is excluded. 

e N, is independent from all other processes N; so that we can factor 
expectations when calculating the risk-neutral discounted payoffs 
for the tranches and the indices. 


Concerning recovery issues, in the dynamic loss model recovery 
can be made a function of the default rate C, or other solutions are 
possible — see the recovery discussion in section 4.1 and also Brigo 
et al. (2007) for further discussion. Here we use the above simple 
methodology to allow losses of the pool to penetrate beyond (1 — R) 
and thus severely affect even the most senior tranches, in line with 
market quotations. 

It is worth presenting an example to further clarify the GPL for- 
mulation. Let us take an example considering a Generalized-Poisson 
Loss (GPL) model where there are three independent Poisson pro- 
cesses N1, N2, N3 multiplied by sizes 1, 10 and 40, with constant de- 
terministic intensities A; = 4%, à2 = 1.5%, 43 = 0.75%. Assume a 
40% deterministic recovery, and take one scenario where the three 
processes each jumped one time before the current time. The loss of 
the standardized DJ-iTraxx pool of 125 names associated with this 
scenario with one jump of N1, N2 and N3, respectively, would be 
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Figure 7.21 Default counting Mr probability distribution associated to a GPL 
featuring three independent Poisson processes with jump amplitude a; = 1, a2 = 
10, a3 = 40 with constant deterministic intensities à; = 4%, à2 = 1.5%, 23 = 
0.75%. 


0.48%, 4.8% and 19.2%, respectively. The expected loss in one year 
for the uncapped loss process is 


0.2352% = 0.0048 0.04 + 0.048 0.015 + 0.192 0.0075 


In Figure 7.21 we plot the resulting default counting distribution 
on four increasing maturities: 3, 5, 7 and 10 years. We note that the 
modes on the right tail of the distribution become more evident as the 
maturity increases. These bumps are the corresponding probabilities 
of a jump of a higher amplitude occurring. 

From the distribution we notice that the probability is not simply 
assigned to the higher amplitude jumps (10 and 40 defaults in our 
example), but also to a number of defaults immediately above: this 
would be the probability of a jump of both high and lower amplitude 
(but higher probability), with size 1 in our case. 

We further notice that some second-order modes start to appear as 
maturity increases. These are the bumps corresponding to multiple 
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jumps in the higher sizes (10 and 40), smaller probability, GPL 
components. 

Let us now assume that we are to price two tranches: 24.8% 
and 10-19.2%. The expected tranche loss of these tranches, and 
ultimately their fair spread, will depend primarily on the probability 
mass lying above the tranche detachment. 

Bearing this interpretation of the modes in mind, we decided to 
choose the GPL amplitudes by fixing independent Poisson’s jump 
amplitudes to the level just above each tranche detachment consid- 
ering a 40% recovery. This led to the values as, ..., @j9 above for 
our as. With an eye to the richness of shapes of the Implied Copula 
historical calibrated distributions, we added four amplitudes (from 1 
to 4), corresponding to a small number of defaults. 

We now present the goodness of fit of the GPL just outlined through 
history. We measure the goodness of fit by calculating, for each date, 
the relative mispricing: 


saz theor_ sa.s.bid 


g Eor z O 


sas mid 
MisrRa ^? = a Btheor  .4,2,aSk 
ae = Si? 4.8.theor bs 5 4-B.ask 
g-emid 0 0 
0 


0 otherwise 
where So :.theor is the tranche theoretical spread as in equation 
(2.4), where in the calculation of the expectations of both numerator 
and denominator we take the loss distribution as resulting from the 
calibrated GPL. 

In Figures 7.22 and 7.23 on pages 120 and 121 we present the 
relative mispricing for all tranches, for all maturities and for both 
indices throughout the sample: from March 2005 to June 2009. 

We note that while the 5-year tranches could be repriced fairly 
well in the current credit crisis for both DJ-iTraxx and CDX, the 
10-year tranches calibrations for both indices have been sensibly 
less precise, following the collapse of Lehman Brothers. Recently, 
with the stabilization of credit markets, we see again fairly precise 
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Figure 7.22 Relative mispricing resulting from the GPL calibration: DJ-iTraxx. 


calibration results (i.e. a relative mispricing in the 2% to 4% range) 


for all tranches, both indices and both maturities. 


To highlight where the problems in calibration come from, we have 


grouped mispricings across three different categories: 


e Instrument: We have grouped all tranches, independently of ma- 
turity and seniority, in one group and have put the remaining cali- 
brated instruments, i.e. the 5-year and 10-year indices, in the other 


group. 
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CDX 5-year — Relative Mispricing 
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Figure 7.23 Relative mispricing resulting from the GPL calibration: CDX. 


e Maturity: We have grouped all tranches and indices in two groups 


according to their maturity: 5 and 10 years. 


e Seniority: We have grouped all tranches (leaving out the indices) 
into three categories depending on the seniority of the tranche in 


the capital structure. 


— Equity: Equity tranche for both the DJ-iTraxx and CDX. 

— Mezzanine: Comprising the two most junior tranches after the 
equity tranche. For the DJ-iTraxx this means the 3—6% and 6- 
9% tranches, whereas for the CDX this means the 3-7% and the 


7-10% tranches. 
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— Senior: Comprising the remaining most senior tranches. For 
the DJ-iTraxx this means the 9-12%, 12-22% and 22-100% 
tranches, whereas for the CDX this means the 10-15%, 15-30% 
and 30-100% tranches. 


From Figures 7.24 and 7.25 on pages 123 and 124 we note that 
the GPL model produced calibrations that resulted in a relative mis- 
pricing that was larger (even though of a relatively small magnitude) 
in the period from June 2005 to October 2006. For both the DJ- 
iTraxx and the CDX, the mispricing could be ascribed to the 10-year 
tranches, in particular the equity and mezzanine. 

We also note that from October 2008 to June 2009 (in other words, 
after the default of Lehman Brothers), the GPL calibration resulted 
again in a non-zero but still fairly contained relative mispricing that, 
in this case, could be ascribed to both the 5- and 10-year tranches 
independently from their position in the capital structure: equity, 
mezzanine or senior. 


iTraxx — Relative Mispricing 


12% 


10% 


8% 


6% 


4% 


2% 


0% 
Mar- Jun- Oct- Feb- Jun- Oct- Feb- Jun- Oct- Feb- Jun- Oct- Feb- Jun- 
0 05 05 06 06 06 07 07 07 08 08 08 09 OY 


— 5y all — 10y all 


iTraxx — Relative Mispricing 
12% 


10% 
8% 
6% 
4% 


2% 


0% 
Mar- Jun- Oct- Feb- Jun- Oct- Feb- Jun- Oct- Feb- Jun- Oct- Feb- Jun- 
05 05 05 O06 06 06 07 07 07 08 08 08 09 09 
— tranche — index 


iTraxx — Relative Mispricing 


14% 4 
12% q 
10% 4 
8% 4 
6% 4 i; 
m Aaa 
4% TA t r tt a 
| fol x A ` uh (2 
2%] Love A TEEPA 
ise ¥ Hae 
0% T * 


Mar- Jun- Oct- Feb- Jun- Oct- Feb- Jun- Oct- Feb- Jun- Oct- Feb- Jun- 
05 05 05 06 06 06 07 07 07 08 08 08 09 09 
— eqty — mezz snr 


Figure 7.24 Relative mispricing resulted from the GPL calibration grouped by: 
maturity (5 and 10 years), instrument type (index and tranches) and seniority (equity, 
mezzanine and senior): DJ-iTraxx. 
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Figure 7.25 Relative mispricing resulted from the GPL calibration grouped by: 
maturity (5 and 10 years), instrument type (index and tranches) and seniority (equity, 
mezzanine and senior): CDX. 
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Final Discussion 


and Conclusions 


8.1 THERE ARE MORE THINGS IN HEAVEN 
AND EARTH, HORATIO... 


We have followed a long path for credit derivatives, and CDOs in 
particular, from the introduction of the Gaussian Copula consistent 
at most with a single tranche at a single maturity, to the introduction 
of arbitrage-free dynamic loss models capable of calibrating all the 
tranches for all the maturities at the same time. 

The critics we presented to the use of the Gaussian Copula and 
of compound and base correlation had all been published before the 
beginning of the credit crisis. Thus, the notion that quantitative ana- 
lysts and academics had no idea of the limits and dangers underlying 
the copula model is simply false. There is even a book titled Credit 
Correlation: Life After Copulas, edited by Lipton and Rennie (2007) 
that is the summary of talks given at a conference of the same name 
in 2006 in London. 

Despite these warnings, the Gaussian Copula model is still used in 
its base correlation formulation, although under some possible ex- 
tensions such as random recovery. The reasons for this are complex. 
First, the difficulty of all the loss models, improving the consistency 
issues, in accounting for single-name data and to allow for single- 
name sensitivities. While the aggregate loss is modelled in order to 
calibrate, satisfactorily, indices and tranches, the model does not see 
the single-name defaults but just the loss dynamics as an aggregate 
object. Therefore, partial hedges with respect to single names are not 
possible. As these are crucial in many situations, the market remains 
with base correlation. 
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Furthermore, even the few models that could try to have single- 
name consistency have not been developed and tested enough to 
become operational on a trading floor. Indeed, a fully operational 
model with realistic run times and numerical stability would be more 
than a prototype with some satisfactory properties that has been run 
in some “off-line” studies. Also, when one model has been coded 
in the libraries of a bank, changing the model implies a long path 
involving a number of issues that have little to do with modelling 
and more to do with IT problems, integration with other systems, 
etc. Therefore, unless a new model appears to be greatly promising 
in all its aspects, there is reluctance to adopt it on the trading floor. 
This inertia, when not excessive, is a good aspect of the modelling 
workflow in institutions, since it gives stability to the modelling 
paradigms being used. On the other hand, if this inertia becomes 
predominant, it will stop all attempts at innovation. A balance is 
required, although the correct balance is difficult to achieve, given 
the number of external (often political) factors affecting the decision- 
making processes of banks and other financial institutions. 

Also, while prior to the crisis a vast number of the quantitative 
analysts were working on credit modelling, after the crisis started 
there was a shift of interest towards more traditional asset classes. 
As the CDO market has become slightly unfashionable, many quants 
have moved to other problems. Counterparty risk valuation is invert- 
ing this tendency to some extent, but the modelling effort in “credit 
correlation” has certainly diminished. 

All these issues concurred in an unfinished research area. To this 
day there is no fully tested operationally satisfying and single-name 
consistent dynamic model capable of a consistent calibration of in- 
dices and CDO tranches across capital structure and maturity. It is 
true that, to some extent, there is more to this problem, but this 
does not mean, however, that the quant community was unaware of 
the limitations of copula models before the crisis, as we abundantly 
documented. 

Where does this leave us? Can we dare express a vision of how the 
future will develop starting from this crisis? 
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8.2 ... THAN ARE DREAMT OF IN YOUR 
PHILOSOPHY 


This brings us into the general picture again, where modelling is just 
one of the many elements that may have contributed to the crisis, but 
certainly not the main one. As we broaden the picture, the number of 
factors that have little or nothing to do with modelling, but have been 
affecting the crisis, becomes apparent. This helps us to put modelling 
into perspective as we close the book. 

There are already several articles that are good references on the 
crisis from a perspective that is relevant for this book, and also 
more generally. We cite a few examples quite different in nature: 
Crouhy et al. (2008); Szego (2009), who presents an analysis of 
the crisis along a path starting with the investors “feeding” on their 
houses; NINJAs! continuing with the frauds, partly documented by 
the FBI’; the problems of securitization and of the originate-to- 
distribute system; the search for yield; the reversal in the real estate 
markets; the absence of supervision; and finally the inadequate or 
wrong capital requirement rules and spillover effects. We also refer to 
El Namaki (2009a, 2009b) for an analysis of managers’ performances 
and quality before and during the crisis. 

We follow, however, Steven Shreve in his Forbes commentary (see 
Shreve, 2008), where he says: 


The models [...] determine the prices of those securities and steer the 
traders who make markets in them. Without this guidance, the banks 
might not have touched them in the first place. To prevent a recurrence 
of financial crises, some call for a return to a simpler time, before 
derivative securities and the quants who analyze them — a time when 
investors bought stocks and bonds and little else. Such complaints 


' Defined as: “An informal term for loans made in the USA to aspiring home-owners with ‘No In- 
come, No Job or Assets’ (NINJA). More formally known as subprime mortgage lending, these loans, 
and the subsequent inability of many borrowers to repay them amid rising interest rates and falling 
house prices, were largely responsible for setting off a global credit crunch in mid-2007,” (The 
British Banker Association, see http://www.bba.org.uk/bba/jsp/polopoly.jsp?d=828&a=10646.) 

2 See, for example, the FBI Mortgage Fraud Report, 2007, www.fbi.gov/publications/ 
fraud/mortgage_fraud07.htm. 
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miss the point. When a bridge collapses, no one demands the abolition 
of civil engineering. One first determines if faulty engineering or 
shoddy construction caused the collapse. If engineering is to blame, 
the solution is better — not less — engineering. Furthermore, it would 
be preposterous to replace the bridge with a slower, less efficient ferry 
rather than to rebuild the bridge and overcome the obstacle. 


This simple analogy is a good description of what, in our opinion, 
is needed in the industry. We need better models. This book has 
followed a path to describe some improvements and a modelling 
evolution, still in need of further improvement. But this is only a 
small part of the picture. To put models more into context, let us 
continue with Professor Shreve: 


Energy trading did not disappear with the demise of Enron, and nei- 
ther will mortgage-backed securities after this fiasco. Put simply, 
the bridge between lenders and borrowers will be rebuilt, because 
we need it. It should, however, be built with better engineering and 
greater simplicity than before. Before the collapse, Carnegie Mellon’s 
alumni in the industry were telling me that the level of complexity in 
the mortgage-backed securities market had exceeded the limitations 
of their models. [...] But in most banks, the quants are not the decision 
makers. When they issue warnings that stand in the way of profits, 
they are quickly brushed aside. Furthermore, in addition to better 
engineering, the bridge must not be built this time with the shoddy 
construction material of no-documentation mortgage applications and 
a network of unscrupulous mortgage originators. 


This reflects two other important points: the first is that the quants 
are indeed not the decision makers. The business drives everything, to 
the point that at times deals are agreed even when there is awareness 
of the inadequacy of models. Again, CPDOs and their ratings come to 
mind as a simple example? where one can even find reports pre-crisis 
warning against the dangers of such instruments and the inadequacy 


3 See the Fitch Ratings report “First Generation CPDO: Case Study on Performance and Rat- 
ings”, published before the crisis in April 2007, stating “[...] Fitch is of the opinion that the past 10 
years by no means marked a high investment grade stress in the range of AAA or AA.” This did 
not stop the CPDO business that was quite strong in 2006 and the beginning of 2007. 


Final Discussion and Conclusions 129 


of the ratings addressing them. But such warnings have largely been 
ignored because the decision makers gave priority to the business. 

The second important point is that, indeed, there have been cases 
where mortgage portfolio deals had to be valued with poor docu- 
mentation. In some cases even the content of the portfolio is not 
clear: as we already mentioned in the introduction, the authors have 
seen offering circulars of a RMBS on a huge portfolio of residential 
mortgages where more than 2% of the properties in the portfolio 
were declared to be of unknown type. What inputs can we give to 
the models if we do not even know the kind of residential property 
that functions as the underlying of the derivative? Also, at times data 
for valuation in mortgage CDOs (RMBSs and CDOs of RMBSs) are 
dubious and can be distorted by fraud, as we have seen above. 

And, as hinted above, we should also include in the picture the 
fragility of the originate-to-distribute system, possible supervisory 
oversight and regulatory errors. The reader can find a lot of discussion 
on such issues too, in relation to the crisis. 

We could continue along these lines, yet what we have discussed 
so far is sufficient to illustrate the complexity of the crisis and the 
limited scope and roles that mathematics and methodology have had 
in it. It should be clear by now that blaming mathematics for the 
subprime fiasco, Lehman, Madoff, bail-outs, unemployment and one 
of the worst recessions in recent history, is missing the point. 

We hope this book has contributed to this understanding through its 
methodological journey and first-hand witnessing of the awareness 
some quants had of the models’ limitations. The steps we — and many 
other researchers — took to address such limitations even before the 
crisis started are a good testimony to the awareness and goodwill of 
several quants and academics. 
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(CDOs) 126 
base correlation 24 
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deterministic recovery 108, 111-15 
DJ-iTraxx indices 44, 65-72, 
109-11, 114 
EETL 64, 68, 70-3 
GPL model 119 
in-crisis 106-15 


large pool homogeneous model 
30-1 

mispricings 67-8, 109-13 

NPV 63-71, 73 

numerical results 71—2 

super-senior tranches 65, 72, 106-8, 
110-12 

surface model 106-15 
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see also Generalized-Poisson Loss 
model 


M 
market quotes 15-19 
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